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ABSTRACT 
Open web steel joists are often used as simply supported 
flexural members in order to support roofs and lightly loaded floors. 
The arrangement of web members permits easy passage of heating ducts 
and other services through the joist. Joists also offer saving in 
weight over comparable members having solid webs. A common chord 


member is a hat shaped section. 


This work is an attempt to predict, by numerical analysis, 
the failure of hat section top chord members of open web steel joists. 
A finite element formulation has been applied to study the behavior of 
these top chord members including non linear geometric and material 
effects. Critical members which failed during an experimental in- 
vestigation are modelled and predicted failure loads compared with the 


experimental results. 


The interaction curves predicted by Galambos and Ketter 
for symmetrical I-Sections and hat shaped sections are used to verify 


the finite element programs. 
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CHAPTER 1 


INTRODUCTION 


1.1 Introductory Remarks 


Open web steel joists are often used as simply- 
Supported flexural members. A common type of joist is shown in 
Figure 1.1. Details of a typical top chord section are shown in 
Figure 1.2. The principal axes, x and y, pass through the centroid 


G. The section is symmetrical about the y axis. 


Matiisen (14) has recently carried out a test program in 
which the top chord members of joists failed due to the secondary 
moment resulting from joint eccentricity. The top chord members of 


these joists were type C sections (26) as shown in Figure 1.2. 


The purpose of this work is to attempt to evaluate the 
maximum load carrying capacity of the chord members of these joists 
by modelling the critical members as inelastic beam-columns. The 
technique employed consists of isolating the critical top chord member, 
attempting to impose realistic boundary conditions upon it, and determining 
the peak load from a Tos eere Teer oh finite element analysis. Epstein 
and Murray (4) developed three different finite element computer 
programs for the analysis of large deformations of beams on the basis 
of an objective non-linear formulation analogous to Budiansky's non-linear 
shell theory. The formulation results in a set of equilibrium equations 


which are solved by the Newton-Raphson procedure combined with an 
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incremental loading technique. The virtual work form of these 
equations has been used to derive a finite element representation 

of the problem, in which cubic polynomials form the basis functions 

for both axial and transverse displacements. It has been demonstrated 
through numerical examples (4) that the resulting equations yield 
computationally efficient solutions to beam problems involving very 
large displacements and rotations. The theory is based on an objective 


measure of strains. 


The three different programs, known as ELAST, PLAST and 
PLAST] deal with linearly elastic, elasto-plastic (unsymmetrical I 
section), and tri-linear stress-strain relations respectively. The 
PLAST] program (4) which deals with the tri-linear stress-strain case, 


includes initial residual stresses for a symmetrical I shape. 


The author has modified program PLAST] such than an 
analysis can be carried out for an unsymmetrical I shape including the 
effect of residual stresses with a tri-linear stress-strain curve. 


This program will be referred to as NLHAT. 


The hat shaped sections have been converted to equivalent 
unsymmetrical I sections, as indicated in Figure 1.3 in order to apply 


the program. 


122 Scope 


Chapter 2 of this thesis deals with a review of various 


techniques available for analysing and computing the capacity of beam- 
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columns. Chapter 3 describes in detail the theory for the programs of 
Epstein and Murray, and the author's modifications to include the effect 
of residual stresses for an unsymmetrical I section. The curves of 
Galambos and Ketter (6) are used to verify the validity of the programs 
PLAST] and NLHAT. Interaction curves for different slenderness ratios 
having symmetric I shape are predicted by the programs and compared 
with those in Reference 6. Hat shaped sections are examined in 

Chapter 4. Chapter 5 consists of a study to simulate the failure of 
the critical top chord member of each of three joists in Matiisen's 
test series (14), under different boundary conditions. A summary of 
the complete work, along with conclusions is the subject matter of 


Chapter 6. 
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CHAPTER 2 
LITERATURE REVIEW 
eS | Elastic Beam-Columns 


Beam-columns are those members in a structure which are 
Subjected to combined axial loads and bending moments. Bending may 
be caused by moments applied to the ends of the member, as shown in 
Figure 2.1, by transverse loads acting directly as shown in Figure 


2.2, or by a combination of these effects. 


The following general basic assumptions are made to solve 
the beam-column problems which involve inelastic bending. 

dl. The cross section of the beam remains plane during 
bending, and hence longitudinal strains are proportional to their 
distances from the neutral surface. 

2. The relationship between stress and strain in any 
longitudinal fiber is given by the uniaxial stress-strain diagram of 
the material. 

he Curvature can be approximated by the second derivative 
of the transverse deflection, neglecting certain higher order terms and 


assuming deflections to be small. 


To obtain the exact solution to a small deflection beam- 
column problem it is necessary to solve the basic differential equation 
of equilibrium. The basic differential equation for the elastic beam 


column problem shown in Figure 2.1 is given as (21) 
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where EI represents the flexural rigidity, P is axial load, Mo denotes 


the bending moment at the left support, y is the deflection at any 


section x, 8B is ratio of end moments and L is the length. Using the 


following notation for simplification 


2 M 
a caer me ae (2.3) 
X 


The general solution of this equation (21) is 
eomeem(e Tox 
yer) cos" kx +.B sin kx -95— [il + r 1 (2.4) 


The constants of integration A and B are now determined from the 


conditions at the ends of beam-column. Since the deflections at the 


ends of the beam column are zero, we conclude that 
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Substituting these values in Equation 2.4 we get 


M F ; 
: = 
you 0 fsin k(L-x) + B sin kx _ ft 2% x) (1-8)] 


P Sin kL 


Now from Equation 2.2 
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P fev p 
lesa ea Pe = 9 |B a (2.7) 
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The moment at the corresponding section is then 


Sin a(l - io +B Sind () 


¢) Sin a (2.8) 


Since the maximum stress will be reached at the section of maximum 
moment, it is necessary to define the location of the critical section. 
This occurs when 


an _ [sin a(] - *) + 8 sina (D1 = 0 (2.9) 


which gives 


cos a(] - fe 
Pies arate (2.10) 
cos on() 
Solving Equation 2.10 for x and substituting Equation 2.10 into Equation 


2.8 gives 
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The interaction relationship to determine the conditions for first 


yield (for the case where no residual stress is present) is therefore 


pt ee = 1.0 (2.12) 
ayy 
M 
sta? SecaX® = 1.0 (2.13) 
Ve ley 
Pp M, 
po tye tT = 1.0 (2.14) 
ay, 
where t= Seca ‘i (2.15) 


This is also known as the secant formula for beam-columns. 
It is an elastic limit solution to a problem where the failure load of 
the column is defined as the load which produces the beginning of 


yielding in the fibers subject to maximum compression. 


2.2 Early Work in Inelastic Behavior 


The use of open web steel joists has increased greatly in 
the recent past. However, very little research has been conducted into 
beam-column problems arising specifically in the top chords of joists. 
This chapter is concerned with a review of current techniques for solving 
beam-column problems and some of the investigations into the compression 


chords of joists. The following summary of early work in inelastic 
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beam-column behavior is taken from Bleich (2). 


Ostenfeld (2) first derived design formulas for axially 
and eccentrically loaded columns. His method was based upon the 
concept that the critical load be defined as the loading which first 
produces extreme fiber stresses equal to the yield strength, in a 
manner similar to that of Section 2.1. Karman (15) was among the first 
to consider the buckling of eccentrically loaded columns as a stability 
problem. He called attention to the high sensitivity of intermediate 
and short columns to eccentricity of the superimposed load, which reduces 


the carrying capacity of straight columns considerably. 


Ros and Bruner (2) assumed the elastic curve of the column 
as a half sine wave and based the computation of the critical load on 
the actual stress-strain diagram. Westergaard and Osgood (18) assumed 
a cosine curve as the deflected center line of an eccentrically loaded 
column, thus simplifying Karman's method (see Section 2.3) without 
impairing the practical accuracy of results. Starting from Karman's exact 
concept, Chwalla (15) based all his computations on one stress-strain 
diagram, adopted as typical for structural steel. His laborious work 
brought insight into the behavior of eccentrically loaded columns and 
his exact results served as a basis for the approximate methods. Jezek 
(24) gave an analytical solution for steel columns based upon a simplified 
stress strain curve consisting of two straight lines and showed that the 
results agree rather well with the values obtained from the real stress- 
strain relation. Shanley (2) showed in 1946 that the double modulus 


theory represented a paradox and this led to the wide acceptance of the 
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reduced (tangent) modulus concept. He said that the column begins to 
bend with increasing axial load and while this load is more than the 
tangent modulus load, i.e., the maximum load for which the column can 


remain straight, it is considerably less than the double modulus load. 


2.3. Karman's Technique 


All of the essential ideas for the computation of in- 
elastic buckling loads are contained in the technique of analysis 


developed by Von-Karman (2). 


Considering the rectangular section shown in Figure 2.3, 
let ey and Eo be the strains in the extreme fibers on the convex and 
concave sides of the beam, and hys ho be the distances from the neutral 
axis to the convex and concave surfaces of the beam, respectively. 
Also let A be the sum of the absolute values of the maximum elongation 


and maximum contraction as shown in Figure 2.4, that is 


=i 


= = Aah ay 
ae Eo oto. 


a (2.10) 


os 


where o is the radius of curvature which may be determined by this 
equation, if A is known. The position of the neutral axis is determined 
by the values of Ey and Eo and is, in general, shifted from its position 
of pure bending by an amount defined by the strain er at the centroid, 


caused by the centrally applied load P as shown in Figure 2.4. 


The force acting at any cross section of the beam can be 


reduced to a compressive force P applied at the centroid of the cross 
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section and a bending couple M. The value of P and M can be calculated 
in each particular case from statics by using the stress-strain diagram 
(Figure 2.4). If y represents the distance from the centroid to any 


fiber of the beam (Figure 2.3), the strain at any point is 
c=e,+¥ (paz) 


Rearranging Equation (2.17) y = p(e-e,) is obtained and hence dy=pde. 


Then the magnitude of the compressive force P is given by 
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ody=- bo f .o de = -—f 
2 we aw 


Dividing P by cross sectional area bh, average compressive stress or is 


obtained as 


a 


2 ae i 
oo Vabhir mn ‘ o de (2.19) 


2 
The integral in this expression represents the area under 
the stress-strain diagram, that is, the shaded area of Figure 2.4. The 
area corresponding to compression is taken negative and the area 
representing tension is positive. From Equation 2.19 the value of Ey 
is calculated corresponding to any assumed value of Ey> provided the 
axial load P is known. Alternatively, both ey and E> may be assumed so 


that the corresponding value of P can be calculated. 


The bending moment is given by the expression 
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The integral in this expression represents the static 
moment of the shaded area of the stress-strain diagram (Figure 2.4), 
with respect to vertical axis A-A. Thus the value of M can be 
calculated for any assumed values of Ey and En: Equation 2.21 can be 
represented in the form 


m= EL (2822) 


where E" = i) 
A 


Oo (Bae) de (29e3) 

By varying Ey and Eo in such a manner that O. remains 
constant, E" is obtained as a function of A = E,-Eo = 2 for any given 
value of Ou. The resulting relation is expressed graphically in Figure 
2.6 which was plotted for a structural steel having the stress-strain 
diagram of Figure 2.5 (21). When these curves are used with Equation 
2.22, the bending moment M can be represented as a function of A for each 
value of O. as shown in Figure 2.7. The shape of the deflection curve 
for an eccentrically loaded bar can be obtained by using the curves of 


Figure 2.7 and applying either an approximate graphical method of 
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integration, as developed by Karman, or a numerical method. 


Consider the steel bar (Figure 2.8) with rectangular cross- 
section, aS shown in Figure 2.3, symmetrical about the point 0, and 
having the stress-strain curve shown in Figure 2.5. Karman's graphical 
integration technique (15) is described as follows. Assuming definite 
value of Ey and Eo for the cross section, corresponding values of P 
and Oo. are calculated from Equations 2.18 and 2.19,Using the curves of 
Figure 2.7, the moment M corresponding to this value of oO. is obtained. 
Thus the bending moment and the compressive force for the cross section 
at the center of the beam are determined, and the distance ota . 
locates the line of action of the compressive force (Figure 2.8(b)). 


Next construct an element 0-1 (Figure 2.8(c)) of the deflection curve 
A 
A 
of the bar. The deflection at the cross-section |] is approximately the 


of a small length a, by using the radius o = ~ calculated for the middle 


2 
same as for a flat circular arc. Thus 6, = oe and the bending moment 


M, = P(5,-5))- For the resulting M, we find the corresponding value 
he 
ma 
better approximation is obtained if we repeat the calculation for the 
oto 
first interval by taking the radius ( 5 I) before going to second interval. 


of A, denoted by Ay» from Figure 2.73; and also calculate yee oA 


Using the new radius, construct the second portion 1-2 of the curve and 
calculate the deflection 64. continuing these calculations to the end 


A of the compressed bar. 


The deflection 6 at the end and the eccentricity e of the 


load P corresponding to the assumed values of Ey and E> are determined. 
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Carrying out calculations for several values of ey and Ens and 
selecting these values in each case so as to make P always the same, 
the deflection 6 is obtained as a function of the eccentricity e. 
Repeating this process for a range of axial loads P and bending 
moments M , Karman constructed equilibrium diagrams such as in Figure 
2.8(d).This method forms the basis for Karman's investigation of beam- 


column compression members. 


2.4 Moment Curvature Relations 


Ketter, Kaminsky and Beedle (11) developed a "point-by- 
point" method for computing the moment curvature relationship for wide 
flange columns bent about either axis and in the presence of axial 
thrust. The necessity of computing M-o curves arose because it was very 
difficult to obtain strains from the given P and M as indicated by 
Equations 2.18 and 2.20. These equations require an iterative solution 
to obtain values of strain, but it is easier to find P and M if strains 


are assumed. 


Since strains are small in comparison to the depth h of 
the section, the curvature within the elastic range as shown in Figure 


2.9(a) is given as 
g: hate Ch Bor tontee (2.24) 


For an elastic perfectly plastic member loaded beyond the 


elastic limit, as shown in Figure 2.9(b), we have 
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eq), = iS Se ees 
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6= Ee (2.26) 


It is thus possible to compute P, M and » for any given stress 
distribution by using equations similar to Equations 2.18 and 2.20. 

As moment is increased under a constant axial force, successive stress 
distributions will be similar to those shown in Figure 2.10. Yielding 
first occurs at the extreme fiber on the compression side of the 
specimen. As the loading is increased, yielding penetrates through 

the web on the compression side. Eventually yielding penetrates through 
the tension sides and the fully plastic condition, typified in Figure 


2.10 (d), may be approached. 


For easy computation of the M-¢ curve certain flexure 
conditions are assumed; a typical set of these is illustrated in 
Figure 2.11. Values of P and M are determined from the stress distribution 
for a number of assumed values of o» for each of a set of assumed values 
of yield penetration as shown in Figure 2.11. These values of P and M 
are plotted in non-dimensional form as 1 VS * and 5 VS s— curves 
as shown in Figure 2.12. The M-¢ curve for a Si Pe fed value of i 
can be determined with the help of Figure 2.12. As an example, the value 
P= 0.2 PY is chosen. A horizontal line is drawn for this 5 cutting 
the various auxiliary P-o curves. At these points of Pa opecee tone 
verticals are projected until they intersect the corresponding M- curves 


in the same figure. Connecting these points give the desired M-¢ curve 


for o> 0.2. This is indicated separately in Figure 2.13. Similar 
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curves may be obtained when the influence of residual stresses 

shown in Figure A.] are considered. The influence of varying the 

axial load over a range of 5 on the M-¢ curve is shown in Figure 2.14. 
The dotted curves in Figure 3.14 are comparable functions considering 


the influence of residual stress. 


2.5 Iterative Numerical Integration Technique 


Ketter (10) presents the result of a continuing investigation 
of the ultimate carrying capacity of beam-columns. He examined two 
different types of loading conditions for wide flange beams and the 
results are presented in the form of interaction curves. The following 
steps were adopted for solving any problems by this technique. 

hes Select a particular value of ~ and P for investigation. 

oe The moment-thrust curvature curve is established for 
the given type of material, cross-section and residual stress distribution, 
for the assumed value of P, as explained in Section 2.4. 

On An end moment M, greater than the initial yield 
moment is assumed. 

4. A reasonable deflected shape is assumed as a first 
trial. The deflections defined by Equation 2.6 are assumed. 


Ds Using P from Step 1, M, from Step 3 and y from Step 


x 
1/2 


Spaced intermediate points on the span. 


4, the bending moment [M = M + Ey is computed at a number of equally 


6. For the value of P selected in Step 1 and the bending 


moments computed in Step 5, corresponding curvatures are read, using either 
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solid or dashed lines from Figure 2.14 depending upon whether or not 
one wants to include residual stress effects, at the points along the 
span. 

Whe Newmark's method of numerical integration (11) is 
then used to determine deflections. Assuming curvatures to be constant 
over a half interval ( x) on either side of each intermediate point 
and considering all slope changes to be concentrated at these intermediate 
points, the angle changes are added to obtain the slope between each 
of the intermediate points. The deflection of each point is then 
estimated by a similar process of adding the slope times the distance 
betieen intermediate points. Then these deflections are corrected to 
satisfy the boundary conditions as explained in Step i of Figure 2.15. 
The final deflection is obtained as the sum of the original estimate 
and the correction. If there is yielding, the deflection will not 
agree with the elastic ones computed in Step 4. 

8. Using the new deflections as a second approximation, 
Steps 5, 6 and 7 are repeated and this set of deflections is compared 
with the second one. This process of successive approximation is 
repeated until convergence to the derived degree of accuracy is obtained. 

9. End slope 85 is computed from me deflection of points 
near the support. If the deflection curve of the member within the 
three end intervals is assumed to be represented by a parabola, the 
end slope can be expressed in terms of known deflections by 
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where 54 and 55 are the deflections at the first and second stations, 
and A is the spacing of the stations. 

10. For the same values of ~ and P, assume a series of 
greater values of M, and repeat Step 4 through 9 for each “4 these. 

11. From the results of the above, plot the yw VS 80° 
One such plot for a particular case of a beam column (when = 0) is 
Shown in Figure 2.16. This curve iS non-linear and eventually reaches 
a maximum, after which it starts to descend. The member would become 
unstable and collapse under the combination of the assigned axial force 
P and M at the peak of the curve. This gives one particular point 
on an interaction curve. 

12. Step 1 through 11 yield one critical combination 
of M and P for one particular slenderness ratio ( ~). Many such points 
need to be determined for a given slenderness ratio to define an inter- 
action curve. 

A series of interaction curves may then be plotted for 
different = values and axial force. Galambos and Ketter (6) plotted 
these curves in non-dimensional form, so as to be suitable for any WF 
section, as shown in Figure 2.17. The effect of Soesnehy Muarsce is 
included. Figure 2.17 has been obtained by solving an 8WF31 section 
with the numerical integration technique described above. 

Mason, Fisher and Winter (13) reported an experimental 
investigation of eccentrically loaded hinged steel columns. Twenty four 
eccentrically loaded and six axially loaded columns were tested. These 
test results were compared with the predicted failure loads of Galambos 


and Ketter, and correlation was shown to be very good. 
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2°6 Interaction Equations 


Mason, Fisher and Winter (13) compared failure loads for 
eccentrically loaded columns to two types of interaction formulas. The 


first was 
ro higi| (2.28) 


where P' is the axial load at failure under concentric compression, 
M' is the bending moment at failure under bending only, 
P is the actual load at failure under combined action, 

and M is the bending moment at failure under combined action. 


The second was 


sr +—— i =] (2.29) 
M'(1 - pT ) 
E 
where Pe is the Euler elastic buckling load and all other terms are 


identical to those of Equation 2.28. Equation 2.28 overestimated the 
test result by 15 to 20% and more and hence it is very unconservative. 
On the other hand Equation 2.29 was safe in all cases in comparison to 
test results. Results obtained from Equation 2.14,which is a secant 
formula ,gives very satisfactory correlation with experimental values 

for ultimate loads for all slenderness ratios. This formula is accurate 
only for bending about an axis parallel to the flanges and provided that 
such bending 7s not accompanied by torsion. Equation 2.29 is accurate 


for bending about either axis when torsion is prevented. The amplification 
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factor (1 - pis included in Equation 2.29 to account for initial 
E 


imperfections. 


Galambos and Ketter (6)present dimensionless interaction 
equations to approximate curves similar to those of Figure 2.17 in the 


following form 
= ] (2°30) 


where A, B and C are empirical coefficients that are functions of = 
and the loading condition, PY 1s the column axial load at the full 
yield condition, M is ultimate bending moment in the absence of 
axial load and Re is ultimate load for a centrally loaded column with 


buckling in the plane of the applied moment. 


Hill and Clark (7) suggested following modified formula 


for both elastic and inelastic range 
s+ —————j—  <] (2°31) 


where us is the maximum applied moment, not including the contribution 
of the axial load interacting with deflections. The possibilities of 
buckling in the weak plane and of lateral torsional buckling are also 
considered in proposing Equation 2.31. Considering unequal end moments 
the following interaction equation is used in modified form in place of 


Equation 2.31 for calculating the critical load. 
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Be cries aes vac tae cnr it Jat Het terse) 
P P 
ae (Mo cr Oe P, 
where Be is the smallest buckling load which the member can sustain 


if bending is absent, (M) pis the critical moment in the absence of 


© 
axial load, and Ch is a factor used for adjusting for the different 
end moment ratios 8. The value of Ch is recommended in CSA Standard 


$16-1969, "Steel Structures for Buildings" (25) for different conditions. 


asdf Investigation of Compression Chords of Joists 


(dy ee Toronto Test Series 


Kennedy and Rowan (9) carried out their work to resolve 
conflicts existing in design requirements for the continuous compression 
chords of open web steel joists, their goal being to establish the 
proper effective length for the compression chord which is a critical 
factor in their design. Thus the testing was aimed at establishing 
the behavior of the compression chord. A total of 8 joists were tested, 
6 of these were loaded at panel points of which 3 failed by in-plane 
buckling of the compression chord, 2 by strong axis buckling of the 
compression chord, due to inadequate lateral support, and one double 
angle joist failed by individual buckling of the angles of the 
compression chord. Two joists were loaded at mid points and they failed 
by in-plane buckling. Deflections, joint rotations and member curvatures 
were recorded to find the inflection points for calculating effective 


length factors. Stub column tests were performed on top chord sections only. 
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From these tests the authors arrived at the following conclusions: 

1. Top chords must be adequately braced. 

Zn Double angles must be suitably joined to ensure 
composite action. 

3. The column strength curve can be closely represented 
by a straight line. 

4, The deflection of a joist can be calculated by simple 
truss theory and can be approximated by a 10% increase in simple beam 
deflections. 

Dp The effective length of a chord member depends upon 
its original shape. If all initial deflections are in the same direction, 
then the effective length factor K = 0.65, if the deflections are not 
in the same direction, then K = 0.9. 

6. When the top chord is uniformly loaded, the drop in 


joist capacity should not be more than 10% if the panel length is 24" 


or less. 
dee Ultimate strength,when bending is present can be 
conservatively predicted by 5 * 7 hee). 
p 
Labaie Kansas Test Series 


McDonald (15) studied analytically and experimentally the 
top chords of Open web steel joists as beam-columns loaded uniformly 
along their length. He did a Karman type numerica| analysis, but 
iterated his integration to find the collapse load which fit his boundary 
conditions. Subsequently, he compared his results with a different inter- 


action equation. A computer program based on the Karman technique (Sec. 
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2.3) was formulated, and used to evaluate the buckling strength of 51 
joists. Twenty four of these joists were tested with one point mid 
Span load, 19 were tested with uniform load applied by Ohmart's air bag 
method and 8 were tested in pairs by uniform dead load. The analytical 
Study was carried out so that results of testing and analytical work 
could be compared. This study was mainly directed towards in-plane 


buckling of the compression chord on the basis of following assumptions: 


l¢ The mechanism of failure is primary instability. 

i The Navier assumption of linear strain under bending 
is valid. 

3% The material possesses equal stress strain properties 


in tension and compression. 

4, The displacement boundary condition for the critical 
top chord member can be determined from the deflected shape of the 
joist by considering inelastic effects from a modified beam analysis. 
Inelastic effects were included by assuming a Variable bending stiffness 


(EI) for the joist along its length. 


Considering the parabolic moment distribution that results 
from a uniform loading of the joist, it was concluded that joists should 
possess their elastic EI in the outer third of their length, and the EI 
over the middle third should be that necessary to match the measured 
mid-point deflection. Thus he decennned joint rotation from the 
gross shape of the loaded joist assuming that the rotation of joints would 
not be affected by the uniform load between panel points. The analytical 


model of a top chord member which has been investigated exhibits: 
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(¢ End conditions which change when the joist load is 
changed and are consistent with the assumed shape which the joist takes 
under load. 

yap An applied load P which rotates to remain tangent to 
the curved form of the joist under load. 
On the basis of his investigations McDonald concluded the following: 

de The analytical study agreed well with the test series 
and both indicated that the design formulas of the American Iron and 
Steel Institute (AISI), Steel Joist Institute (SJI), and the American 
Institute of Steel Construction (AISC) will give adequate factors of 
safety. However, the analytical study indicated that the factor of 
safety might be reduced if the ratio of uniform load to axial design 
load in the top chord was increased. 

2% Design techniques which use only the slenderness 


ratio and axial load to proportion compression members are not rational. 


=e The interaction formula he proposed is 
Axial Stress Bending Stress _ 9 go 


Tangent Modulus Stress ch Yield Stress 


4, Von Karman's method of analysis for beam columns 
may be applied to the compression chord of open web steel joists but 
in analytical studies an iterative process should be used to determine 
the indeterminate end moments which produce a deflected shape consistent 


with the boundary conditions. 
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2.8 Large Deformation Finite Element Solution 


Epstein and Murray (3) presented an objective non-linear 
formulation of elastic in-plane beam-column problems and a corresponding 
set of finite element equations. A large deformation theory for in-plane 
beam problems was formulated on the basis of the Budiansky's non-linear 
shell theory. Cubic interpolating functions were used to represent 
finite element displacements. These displacements were introduced 


into-a virtual work form of the equations to yield numerical solutions. 


A computer program was developed to show the applicability 
| of the theory. Programs have also been developed to solve the post 
buckling behavior of beam columns. The solutions for beam-column 
problems carried out herein have been obtained by this technique and, 
therefore, a detailed discussion of these programs is aiven in Chapter 
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CHAPTER 3 


AN INELASTIC LARGE-DEFLECTION FINITE ELEMENT TECHNIQUE 
FOR BEAM-COLUMN TYPE PROBLEMS 


onl Background to Finite Element Formulation 


Epstein and Murray (3) have developed a finite element 
formulation in order to obtain numerical solutions for non-linear 
inelastic in-plane beam-column problems. Reference (4) contains three 


programs and the theory on which these programs are based. 


The methods of Chapter 2 are based upon the approximations 
required to derive a second order differential equation. However, 
Reference (4) uses an exact expression for the curvature. This theory 
has been developed to solve problems of very large deflections in 
the post-buckling region and in this sense is more elaborate than 
required to determine maximum load carrying capacity. Use is made of a 
Newton Raphson procedure combined with an incremental loading technique 
on a set of total Lagrangian equations in solving the problems. Since 
a large deformation theory should include large rigid body motions 
as a special case, products of small terms have been retained so that 
an objective set of equations is obtained. This chapter contains a 
summary of Epstein's formulation (3) and the resulting programs. The 
basic assumptions made are: 

1. The x,z plane (Figure 3.1) is a plane of symmetry 

of the beam and of the external loads. (The z axis 


contains the axis of the beam before deformation. ) 
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ae Normal sections remain plane, undistorted and 


normal to the axis of the beam after deformation. 


Referring to Figure 3.], the position vector of a point on 


the axis of the beam before deformation is given by 

yr = 2k (ore) 
and for the same point after deformation 

r Eee ri | (222) 
in which displacement vector V can be written in component form as 

v = Uj + wk (323) 


where i and k represent unit vectors in the coordinate directions. 


The natural base vector along the deformed axis iS given by 


WwW > 
i 
Qf 
NI > 


= (1 ¢+w')k tu'i C34) 


where symbol ()' and (*) are used to denote z-derivatives and 
magnitude after deformation, respectively. The unit normal vector 


after deformation is 


A Se le (3.5) 
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The Green strain measure of the axis is given by 
(3.6) 


Substituting Equation 3.6 into Equation 3.5, N is expressed as 


-u'k + (1 + w')i (3.7) 
2et] 


Considering displacements of points not originally on the axis of 
beam in order to construct the two dimensional strain tensor (the 
third dimension, y being irrelevant), and using N to denote the unit 
normal vector before deformation (N=i), the position vector of a point 


in the plane of symmetry x,z is given before deformation by 


R=r+xN (32.8)) 


In the view of assumption 2, the same point is identified after 


deformation by the position vector 
Reel och x’ N (3.9) 


The components of Green's strain tensor are given by the matrix equation 
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[el = 5 : (3.10) 


The vector derivatives appearing in Equation 3.10 are 


evaluated as 


dR *: 

Agha N (Ser) 
Guede dn ay di 

ret dale wets a adzZ (3.12) 


where g is given by Equation 3.4 and N is a vector tangent to 


the axis after deformation and expressible, therefore, as 


gy ein (3.13) 


The scalar » is a measure of the curvature of the deformed 
axis and it identifies with the “physical” curvature (i.e. || = Fae 


where s measures length along the curve). 


Substituting Equation 3.13 into Equation 3.12 


ae esta (3.14) 
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Using Equations 3.11 and 3.14 yields Green's strain tensor as 


defined in Equation 3.10 


Le] = (3.15) 


Since all but the first component of the strain tensor 
vanish, as a direct consequence of the above hypotheses, e€ is used to 


denote non vanishing component of [e]. 


3.2 Variational Equilibrium Equations 


The internal virtual work (IVW) is obtained by integrating 
the following integral extended over the volume of the beam after 


deformation. 


“~ 


IWW=f fo6edA de (3.16) 
Q A 


Here o is the stress component corresponding to ce, & is 
the lengtn of the axis after deformation and A is the area of the cross 


section of the beam. 


Since e and $ are independent, the variation Se is obtained 


from Equation 3.15 as 


6 e@ - 2(1-ox)x (e + =) S¢ (3.17) 
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Defining a modified curvature as 
@ = (2e + 1) ¢ (3.18) 


it can be seen. that for smal] strains reduces to 9. Substituting 


the variation 6¢ from Equation 3.18 into Equation 3.17, 
6 € = ((1-6x) + ox(1-ox)) 6 e - x(1-ox) 66 (3.19) 


Substituting Equation 3.19 into Equation 3.16 and defining stress 


resultants as 


N* = f o(1-ox)dA (e720) 
A 
M* = - f o(1-ox) xdA (3.21) 
A 
Equation 3.16 can be written as 
IVW = s {N*¥Se + M*(do - gdde)} de (3722) 
Wi 


Introduce modified variables, defined as 
N= 2e+] N* (35 23) 
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into Equation 3.22 to avoid the following two inconveniences, 
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es Integration extends over the deformed axis, and 
(ag: $ and are not rational functions of displacements. 
AV) AV) 


Again, for smal] strains N and ¢ reduce to N* and o 


respectively. Define another variable 


Steele aw oat | (3.25) 
where 
AV) 
gx = 2 = J 2e+1 9 (3.26) 
2et] 


Introducing these modified variables expressed by Equation 3.23 to 
Equation 3.25, Equation 3.22 can be written as 

yy 
N 


AV) 
IVW = Ss (N de + M* 6o)dz (3.27) 


Me 


We can note the following about Equation 3.27: 


1. It is exact, in the sense that it contains no 
approximations whatsoever concerning the magnitude 
of displacements or strains. It is equivalent to 
the exact three dimensional internal virtual work 
expression, provided the two basic assumptions of 


the beam theory are valid. 
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2 Integration extends over the original configuration 
of the beam axis. 

SY All the strain variables appearing explicitly in 
Equation 3.26 are rational expressions in the dis- 
placements u and w and their derivatives. 


Two types of loads are considered: 


1. A force per unit undeformed length of axis. 
f = t. Kock 3 i (3.28) 
lbp A force per unit deformed Jength and acting 


normally to the axis after deformation 


cay 


P=pN (3.29) 


where P is an external force and p is its magnitude. The total 


force per unit undeformed length is then 
q = /(2e+1) PN + fk + ue (3.:30)) 
or, according to Equations 3.5 and 3.7, Equation 3.30 can be written as 


q = (f.-pu')k + (f 4p(14w'))i (3731) 
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The external virtual work (EVW), corresponding to the two types of 


external loads discussed above, is given by 
EVW = f q 6 vdz 
L 
Suh, {(F,-pu' ) dw +(# +p(14w') )dutdz (3432) 
y 


The total variational equilibrium equation is obtained upon excluding 
the end forces, by equating the right hand sides of Equations 3.27 and 


3.32, so that 


AV, 
% 


AV) 
SIN de + M*Sp - (f,-pu')dw-(F,-p(1+w') )Su}dz=0 
Q 


(3733) 
Substituting 
e=w' + wi? + ut?) (3.34) 
and 
& = u"(14w') - utw" (3.35) 


into Equation 3.33 and integrating by parts, we obtain 
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the Euler equations corresponding to Equation 3.33 as well as the 
natural boundary conditions. The resulting differential equations of 
equilibrium are 


AV) 


[N (l4w')]' + [M*u"]' + [M*u']" + f.-pu' = 0 (3.36) 
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3.3 Discretization 


A set of algebraic equations is obtained by dividing 
the beam axis into elements and allowing the displacements in each 


element to belong to a finite parameter family of functions. 


Third degree volynomials are used to define the displacement 
functions for the given degrees of freedom as shown in Figure 3.2. The 


basis functions are expressed as 


fie mo 2B 1BGo xt (3.38) 
ioe (ce cuper any eh hy (3.39) 
fy = - 20° + 30° (3.40) 
fa = (ce hy (3.41) 
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in which ¢ is a non-dimensional] coordinate which varies from zero 

at the left end to 1 at the right end of each element, and where 

hy is the length of the element number I. The coefficients of these 
functions in the linear combination are the displacements and their 
first derivatives at the ends of the element. The nodal displacements 


Uy and Ws and their derivatives Ur and We are shown in Figure 3.3. 


There are four degrees of freedom at each node in the 
element denoted by U, W, U' and W'. Making use of Equations 3.38 
to 3.41, the following expressions are obtained for the displacements 


h 


within the {t element and for their first and second derivatives in 


terms of the notation shown in Figure 3.3 
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= Uocperyag ll * Mocreayek? ky? 57 Yocrery4g Mecr-1)+k 


atte (3.44) 
Se = Wo ryyagh g*CWo( 7-1) 4jWacr-1)+k * Y2(-1) 43 
Np rales 
oeoti-1ytk) 3. k (3.45) 
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and the internal virtual work is evaluated as 
2% N 

IVW =f (N de + M*6d)dz = = 
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* Wo (reayek Vacreryg * Yacreryeg Me(r-1)44) 
54 (Maer 2 Maat 3 (3.47) 


in which the following notation is introduced, 


Gt a 
NAA Ee IN fs cde (3.48) 
ani j 
Q 
N¥**, = f iy Lp bea temas a (3.49) 
jk 5 a omeek 
Le 
M¥*. = Sf M*f" dz (3.50) 
dae j 
6 
KKK = Mee" f! ; 
M 3k J 3 j mls (3251) 


dependence on I being understood. 
The external virtual work (assuming dead load forces and moments are 


concentrated at the nodes only), is given by 


N 
EVW = A X Uy yyy * MyL(1 + W545) 8Up s40 
~ Yo ae2 Saga] ee 


in which X 3 zy and My are transverse force, axial force and moment 


at node J. 


The equilibrium equations associated with Equations 3.48 to 
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3.52 are obtained by assuming sequentially that all but one of the 
variations vanish. Thus, for instance, assume that the only non- 
vanishing variation is Sup = 5U5) 47° In Equation 3.47 all the indices 
affecting SU are of the form 2(I-1)+j, the only possibility for a term 


containing 8U5 (1-1) +j not to vanish is that 

2(I-1)+j = 2L + 1 (3.53) 
or, solving for 1, 

r=L- 33 (3.54) 


Since j is a positive integer not greater than 4, it must necessarily 
be either 1] or 3, resulting respectively in I = L +] and I=L. There- 
fore, for each non-vanishing variation considered, there are in general 
only two surviving contributors to the overall summation of Equation 
3.47, namely those corresponding to the elements adjacent to the node 
considered. Only one element contributes to the summation when first 


and last nodes are considered. The equilibrium equations are 


+ N*** M** 


2 
eo 5(2k-1) Yo(a-Ker)4g * Moca-Ker) +g" (2K-1) 5 


4 Mmm oT) = Xj (3.55) 
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2 
kk ek 
Bae kel aK * NS (oK) Yocaker 4g * Mo(a-ker +5 ** (2K) 
“MEHH ony) ] = Mg(1 + Woyyo) (3.56) 
2 : 
kk KKK 
Bhs Ke2 EN okey tN (oK1)Mo(a-ker 4g * Y2(a-ner) eg M5 (24-1) 
- M*xx* = j | Z 
(2k-1)5)4 Z 63757) 
2 
* a 
ya Keak *oK tN*** 5 (Ky Wocg-ker) * Yocg—Ker)4g M5 (2K) ~ M* (2K) 5) 
AST gaade2 Bab) 


where J = 0,1,2,...N. 


A boundary condition of displacement results in the 
impossibility of assuming a non-vanishing variation for the particular 
displacement prescribed. The corresponding equilibrium equation must 
therefore be dropped and substituted by the specified value of the 
unknown. These conditions do not necessarily have to be imposed at 
the ends. They can be specified at any nodal point, allowing the 


solution of a continuous beam. 


Due to the particular choice of stress and strain measures 
in the analytical formulation, the resulting discrete equilibrium 
equations (for the case of a linearly elastic material) are cubic poly- 
nomials in the nodal displacements and their first z-derivatives. 

The Newton-Raphson procedure is adopted to solve the set of equations 


(Equations 3.55 to 3.58). The Jacobian matrix is calculated in which 
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the elements are derivatives of the equilibrium equations with respect 
to the unknowns. This matrix results in a band matrix of width 12. 

It is partitioned into submatrices of 4x4 which are arranged in a tri- 
diagonal manner. Potter's technique (20) is used at each step to solve 


the set of linear equations. 


3.4 Computer Programs 


Three different computer programs have been developed on 
the basis of the theory explained above. They are referred to as ELAST, 
PLAST and PLAST1 and their scope is as follows: 
(a) ELAST deals with the linearly elastic case. 
(b) PLAST deals with an ideally elasto-plastic-material. 
Plasticity effects as such, however, are not 
accounted for. The cross sectional shape may be 
a combination of three arbitrary rectangles, that 
is, an unsymmetrical I. 
(c) PLAST] deals with a tri-linear ‘elastic' material. 
Initial (residual) stresses are included but the shape 


has to be a symmetric I. 


Double precision is used and all the programs include an 
automatic recording of the present step and the last converged solution 
for added flexibility in restarting the iterative process. A more 


detailed description follows. 
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eel Wa Program ELAST 


The program ELAST has the following limitations. 
» A constant element length is assumed. 
2. Constraints on degrees of freedom may be imposed 
at any node, but elastic supports are ruled out. 
3: Loads (forces and moments) are assumed concentrated 


at the nodes. 


For details of data preparation, Reference (4) should be 
consulted. Two types of output are generated by the program, viz., 
printed and permanent disk files. The printed output consists of 
intermediate results (with no explanatory headings) and final results 
(displacements, forces and co-ordinates). The output on disk is 
intended to allow for the restarting of the iteration process. The 
last trial vector is always stored in file 15 (erasing the previous 
vector). It is also convenient to have at hand the last converged 
solution. This may be accomplished by specifying NFl=1 (see a descrip- 
tion of data preparation in Appendix D) which will cause the program 


to write the solution on file 16 only if convergence is achieved. 


3742 Program PLAST 


The equations of equilibrium are the same as for ELAST, 
but yielding effects imply the existence of an upper bound. For the 
external load, an augmented structure technique (23) on this load 


carrying capacity is used. For post-buckling behavior, the augmentation 
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has been realized only by means of linear horizontal, vertical and 
angular springs at the last node (right end) of the beam. The 
limitations on other procedures, including output are the same as 
for ELAST except that only the degree of freedom associated with 


the last node can be elastically supported. 


S453 Program PLAST 


PLAST] is a generalisation of PLAST in two respects: 

(a) It considers a trilinear constitutive equation, 
and 

(b) It allows for the inclusion of initial stresses 
(linearly distributed along the flanges and 


constant on the web). 


It is assumed that the walls are thin enough so that 
stresses can be considered constant through the thickness and the 
integration of stresses over the cross section is not complicated. 
Analysis of symmetric I sections only is allowed. The limitations are 
the same as in PLAST but final results in the output will include 
positions where changes in modulus of elasticity occur (in several 


sections along each element). 


The integrals appearing in the equilibrium equations are 
evaluated numerically with the use of Simpson's rule over the length 
of element. Elements of the Jacobian matrix are evaluated by numerical 


differentiation. 
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3.5 Modifications to Programs 


To achieve a generalisation to other shapes with 


consideration of a trilinear constitutive equation and inclusion of 


initial stresses, the program PLAST] has been modified by the author. 


The flow chart is essentially the same as that in Epstein and Murray's 


report (4). 


8.521 


Modification has been made to the following subroutines. 


lhe 
2. 
3 


MAIN 
ABCG 
FORMOM 


Modification to MAIN 


MAIN has been modified such that: 


(a) 


(b) 


Combinations of linear horizontal, vertical and 
angular springs can be used at any node in the 
entire beam. 

A residual (initial) stress pattern for unsymmetric 
I section has been calculated from the equilibrium 
considerations. The stress on the top and bottom 
of the web and in the narrow flange are determined 
from the given initial compressive stress at the top 
of the wide flange. Initial stresses are assumed 
constant on the narrow flange and linearly varying 
along web and the wide flange as shown in Figure A.2. 


The detailed calculation of the stress distribution 
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is indicated in Appendix A. This pattern of 
residual stresses is intended to simulate those 


arising in a hat section. 


Tne centre of gravity of an equivalent unsymmetrical I- 
section used to simulate a hat shaped section, is also calculated by 
the program. The program reads all] data; including the initial trial 
vector of displacement which is given either on file 15 or file 16. 
Output has been modified so that the program prints all initial and 


final yielding values for both of the flanges and the web. 


3.5.2 -- Modifications to ABCG 
j 


Subroutine ABCG evaluates the non-vanishing band (width = 
12) of the Saboran matrix, by means of DNMP, as well as the right 
hand side, by means of NMP, for use in the Newton-Raphson procedure. 
Thus subroutine ABCG formulates equilibrium equations and its 
derivatives. It has been modified such that augmenting springs can 
be applied at any nodes to get the descending part of load-deflection 
curve without altering the band character of Jacobian matrix. The 
band of the Jacobian matrix is subdivided into a tridiagonal partition 
of 4x4 submatrices for use in POTRS. Force and moment vectors are 


also modified since spring stiffnesses have been applied. 


CaS: Modifications to FORMOM 


Subroutine FORMOM has the function of calculating the axial 


force and bending moment stress resultants from the given two measures 
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of strain ( e and ¢ ) for the given symmetric I-section and for a 
given trilinear constitutive equation. It is modified such that it 
calculates the axial force and bending moment stress resultants for 


a given equivalent unsymmetric I-section. 


The stress and strain diagram and the system of co- 
ordinates are exactly the same as indicated in Reference 4. The 
technique of computing distances to the initiation of yielding in 
the web (y values) has been changed because the residual stress is no 
longer constant in the web but varies from compression at the top 
to tension at bottom. Derivations are shown in Appendix B. With 
these modified values of y, values of axial force and moment stress 
resultant are computed. Once this subroutine was tested and modified 
it was attached to the main program. MAIN with all its subroutines 
is stored in object form in a permanent file called NLHAT. The function 


of each subroutine is described in Reference 4. 
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CHAPTER 4 


ILLUSTRATIVE APPLICATIONS OF PROGRAMS 


4.1 Testing of Program ELAST 


In order to verify the capabilities of the programs 
developed by Epstein and Murray (4) to predict the behavior of beam- 
columns, the programs have been used to solve problems whose solutions 


are already known. 


The first problem chosen to verify the programs is 
the elastica (21), shown in Figure 4.1, which involves the elastic 
post-buckling behavior of a simply supported beam-column. A small 
transverse load was added in addition to the axial compressive force 
P to act as an initial imperfection. Results so obtained are tabulated 
in Table 4.1, where they are compared with the elastica solution (22). 
The plot of Figure 4.1 shows good agreement between these two sets 
of results in the post-buckling region. The number of elements used 
was four for one half of the column taking account of symmetry. The 


method of data preparation is described in Reference 4. 


The solution to this problem by program ELAST was obtained 
by Epstein and Murray in Reference 4. It is repeated here to demonstrate 
the capability of the program to handle very large geometric non- 
linearities. Since the geometric formulation is identical for al] 
programs contained in Reference 4, the capabilities of the programs 
PLAST and PLAST! to represent similar types of deformation have been 
established. 
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4.2 Testing of Program PLAST 


Use of program PLAST is made herein to analyse for plastic 
post buckling behavior. The solution for a W-8x3] simply augeonben 
beam of 210" span, subjected to the action of a constant transverse 
load Q at mid span, and to a varying uniform compressive force P, is 
Shown in Figure 4.2. A sketch of the model, taking account of 
symmetry, iS also shown in this figure. The points plotted are 
tabulated in Table 4.2. The transverse load Q was taken as 3413 lbs, 
so as to Simulate an initial imperfection of 0.2". The boundary 
conditions at the left end are u=w=o, and at the right end u'=o. u 
and w are the vertical and horizontal deflections in the beam, respectively 


and u' is the slope at the end of the beam. 


The half beam was divided into four finite elements. 
The yield strain was assumed to be 0.12% and the modulus of elasticity 


to be 30 x10° psi. 


For the first trial an initial vector of zero is assumed 
with a compressive load of 200,000 lbs. It can be seen from Figure 
4.2 that this load is quite close to the carrying capacity of the member. 
Since no plastic penetration is present, convergence was achieved very 
rapidly. Since yielding effects imply the existence of an upper bound 
for the external load, an augmented structure technique (23) is used. 
Augmentation is realized by means of an application of springs in the 
beam-column. Hence for the second load step a horizontal spring with 
a stiffness of 5 x 10° lb/in was added, which corresponds to roughly 


double the stiffness of the beam in that direction. The spring stiffness 


ae oan eee a 
oH "ae Se DS 
ef A is m : 


‘ NY ‘ ae 


: Pear" & <, 
nae ieee - 


— 


biosiee: J aie 
a ey yy i ae yeni aE 7 


B com'y tog aig oe ta.aeh, ee Stan br 
Viavisaaqesy Mead edd) mt et ° | Sis Lesley sit ve | | 
iiicitie Gh | df Yo Oe sate saoteéat a ‘whe | 1 7 
st i i ak ae 
enters Zufabom oft bie ae See ws sy Us | 
ae ican ore ea e i ¥ By ipa hectigee 


ni a > Lok 7 Es 


ie, Fee apie he 


nia bebo | 


' <3) 
alee 


a Ea 


me i rian ee 
frpssterteihagette 
analy yaa 7 . 


48 


is doubled to achieve more rapid convergence and a smooth load 

deflection curve. This means that the total stiffness is approx- 

imately three times the stiffness of the member and a force of 600,000 lbs 
should be applied to maintain the augmented structure in its equilibrium 
configuration. As the load is increased, more and more of the 

increment goes into the spring until] net force in the original structure 
diminishes. The total load and net load are shown in Table 4.2 as 

well as in Figure 4.2. If a very high load is adopted, the computer 

time limit may be exceeded, with the intermediate results showing a 
divergent behavior. This indicates that the load step is too large, 


and a smaller step has to be adopted. 


Since the contents of file 15 represent the last trial 
vector, NF = 2 is specified for the subsequent run. This will cause 
the program to take the last converged solution as the first trial 
vector. Convergence may then be achieved at a smaller load step. The 
logical thing in using this program is to give a very small increase 


in load steps when plastic effects are about to begin in the section. 


The computer time for convergence is relatively high when 
the load is either at the maximum of the load deflection curve or 
when the structure has passed its maximum carrying capacity. This occurs 
because it requires a greater number of cycles for convergence. The 
plastic penetration is very deep and iS a maximum at the last equilibrium 
position of this example.affecting more than 70% of the beam depth as 


shown in Table 4.2. 
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This example problem for program PLAST is that used 
by Epstein in Reference 4. It illustrates the technique of applying 
Springs to augment the stiffness of the structure in the descending 
region of the load deflection curve, in order to maintain a positive 
stiffness of the augmented structure. The presence of initial 
stresses smoothes the peak of the moment rotation curve. Also, 
the strain-hardening phenomenon has its effect on slowing down deform- 


ations past the maximum load. 


4.3 Testing of Program PLAST] 


The end rotation for the solution of a beam-column problem 
with constant axial load and increasing end moment is tabulated in 
Table 4.3 and plotted in Figure 4.3. The residual stress patterns 
that of Figure A.1 with Onc of 9900 psi (-30y). The geometric and 
loading conditions are illustrated in the sketch in Figure 4.3. In 
this problem Oy = 33 ksi and it is assumed that at a strain of 0.24% 
the strain hardening process begins with a modulus of 6 x 10° psi. 
Ordinates of points where changes in elasticity occur are also 
obtained in the computer output. The maximum point on this curve 
corresponds to the highest end moment which this column can support, 
thus giving the point on the interaction curve of Figure 4.4 with 
>= 0.8 and = = 40. Twelve such points, that define the maximum 
eae for various values of ~ and 5~ have been obtained and ane. 


plotted in Figure 4.4. Points so obtained have been compared with the 


curves of Galambos and Ketter (6), shown by the solid lines in Figure 
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4.4, which include the effect of residua] stresses. The points 

used for Galambos and Ketter's curves are tabulated in Table 4.4 and 

a numerical comparison is given in Table 4.5. These results differ by 
from 0 to 8% from those of Galambos and Ketter and indicate that 


PLAST] gives reliable results for beam columns. 


4.4 Applications of NLHAT 


PLAST] was modified, as explained in Section 3.5, to 
obtain the program NLHAT. It is difficult to obtain solutions in 
the literature for the verification of NLHAT since little attention 
has been devoted to the inelastic instability analysis of unsymmetric 
sections. However, the modifications are essentially concerned with 
the pattern of residual stress and the support conditions and therefore 
the fundamental capabilities of PLAST] have been retained. The ability 
of NLHAT to predict the behavior of an unsymmetric section can be 
verified by testing for the plastification of the cross-section under 
the action of combined axial load and moment, since the conditions for 
complete plastification can be predicted by hand computation. Therefore, 
moment curvature curves were obtained from the program for a very short 
segment of the section and the fully plastic conditions from these 


results were compared with those obtained by hand computation. 


The moment-curvature response for a segment of aC type 
chord section, with Ls = 33000 psi and o = 9900 psi, is shown in 


Figure 4.5 for various ratios of  - In order to simulate this section, 
y, 
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the unsymmetric I - section shown in Figure 1.3 was used. Section 
properties of the simulated section and the chord section are given 

in Table 1.1. In deriving the dimensions of the simulated section 
the objective was to approximate the area and the moment of inertia 

as closely as possible. The length of the segment used in the 
computer analysis was 2 inches and therefore the plots in Figure 4.5 
actually apply to a hat section column with => 4.8. Note that two 
moment curvature curves are obtained for each fixed value of P, one 
for moment subjecting the wide flange to compressive stress (Table 
4.6) and the other for moment subjecting the narrow flange to 
compressive stress (Table 4.7). The fully plastic moments, for the 
various values of 5 obtained from the maximum moments on the 
moment-curvature Peres of Figure 4.5, are tabulated (for ~. =0) in 
Table 4.8 and plotted as the two solid interaction curves in Figure 
4.6. (Moments producing compression on the narrow flange are considered 
positive and those producing compression on the wide flange are 
considered negative.) Hand computations for the fully plastic inter- 
action curves of this section are carried out in Appendix C and plotted 
as the two dashed interaction curves in Figure 4.6. A comparison of 
these curves indicates that the finite element program underestimates 
the stress resultants required to produce fully plastic conditions for 
both signs of moment, the maximum error in Y- being about 10%. This 
indicates that instability loads predicted y the finite element 


technique may be considered conservative. 


For another comparison of the fully plastic interaction 
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curves, the curve for a symmetric 8WF3] section, as plotted by 
Galambos and Ketter (6), is also shown in Figure 4.6. It should be 
noted that both the finite element and hand computed interaction 
curves for the hat section straddle that for the symmetric section. 
Since Galambos and Ketter's curves (6) are based on the behavior 

of an 8WF31, this indicates that their curves are not conservative 
for unsymmetric sections, if the applied loads produce compression 
on the narrow flange. This is demonstrated below for two different 


Q : 
— ratios. 
fe 


Stability analysis, as described in Section 4.3, were 
carried out for the simulated hat section for = ratios of 40 and 120. 
The results of these analyses are tabulated in Table 4.8 and plotted 


as the solid lines in Figure 4.7. 


Galambos and Ketter's curves (6) are also plotted in 
Figure 4.7 for the same slenderness ratios. These results indicate 
that the moment capacity of this hat section may be over estimated 
iv 


by as much as 35%, for fo 40, by the curves based on symmetric sections. 


Interaction curves for unsymmetric sections may be 
conveniently displayed as shown in Figure 4.8, where the lack of 
symmetry due to the sign of the moment is apparent. The values plotted 
in Figure 4.8 come from a series of analyses which are tabulated in 


Tables 4.9 and 4.10. 
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CHAPTER 5 


A STUDY OF TOP CHORD BUCKLING OF JOISTS 


bel Description of Experimental Project 


Matiisen (14) has investigated the effect of joist 
eccentricity on the behavior of joints by performing tests on a 
series of joists. He derived the following conclusions: 

Tes The primary effect of joint eccentricity is to 

induce moment in the chord. 

ria These moments, in the elastic range, can be 

predicted with good accuracy by modelling the 


joist as a frame. 


In his test series five of the seven joists failed due 
to top chord instability. Two of these five had a type of joint 
which transferred little moment into the chord. The remaining three 
had a type of rigid welded connection which was capable of transferring 
substantial moments into the chords. In these three joists, designated 
as AXO1, AX02 and AX05, failure was induced in a top chord member 
which did not have the maximum axial force. Thus, it is apparent that 
the moments in the critical members had a significant influence on 
the benavior of these members and contributed to their collapse. This 
chapter is concerned with an attempt to model the critical member in 
these joists as a beam-column, to predict the failure load of this 
member on the basis of an analysis of the model, and to compare the 


predicted capacity with the failure load observed in the tests. 
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The top chord member of joists AX01], AX02 and AX05, which 
failed due to inelastic instability, is shown as member 3T in Figure 
5.1. The measured dimensions of these joists are shown in Figures 5.2 
and 5.3. Concentrated loads were placed between top chord panel 
points to simulate uniform loading in joists AX01 and AX02, whereas 
in joist AX05 loads were applied at the centre of each of the top chord 
joints. In addition to testing the joists, material testing and elastic 
plane frame stress analyses, which computed joint rotations, displace- 
ments, axial forces and moments at the ends of each member, were 
carried out by Matiisen (14). The pertinent results which effect member 
3T are tabulated in Tables 5.1 and 5.2. The member and joint numbering 


corresponds to that in Figure 5.1. 


In joists AXO1 and AX02 there are large moments generated 
in member 3T due to top chord loading, chord continuity and joint 
eccentricity. However, the moments in member 3T of joist AX05 are due 
primarily to joint eccentricity and continuity effects. A summary 
of the failure conditions for these joists is given in Table 5.3. All 
joists had C-section top chords. 

The present study attempts to predict the load at which 
instability failure occurs by modelling the critical members as 
beam-columns with approximate boundary conditions, and compares 


predicted failure loads with those shown in Table 5.3. 


5.2 Description of Models and Boundary Conditions 
Member 3T was, at various times, analysed with three 
different kinds of boundary conditions, namely, force, spring and 


displacement boundary conditions. AX05 was analysed under all three 
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boundary conditions, wnile AX0] was analysed under spring and force 
boundary conditions, and AX02 was analysed under only the force boundary 
condition. Force boundary conditions were also applied to models 

that included the joint as a part of the member, for all the joists. 
The program used for the analysis was NLHAT, described in Chapter 3. 

The method of data preparation and a listing of the program are given 


in Appendix D. 


A typical finite element mode] of a top chord member 
with n elements, as idealized herein, is shown in Figure 5.4. The 
length of the member between panel point joints is denoted by &. If 
a joint is included in the model, the length of this joint is denoted 
by hi. Loads on the model may consist of the axial force (shown as P), 
loads applied at the panel point joint (shown as Xj, Zj and Mj), and 


sets of interior loads (shown as Xi, Zi and Mi). 


A summary of the geometric dimensions, the joint loads, 
and the interior loads, which are applicable to the six models used 
herein, iS given in Table 5.4. Unfortunately, the program required a 
model with elements of equal length. The number of elements was chosen 
to obtain the joint loads as close as possible to a nodal point. For 
those models which include the joist joint, a J has been appended to 
the joist designation. For interior loads, the loads were replaced 
by statically equivalent forces at the two joints which straddled the 
actual point of application of the load on the joist. These statically 
equivalent Tdads are also Shown in Table 5.4. The technique of obtaining 


joint loads will be more fully discussed in subsequent sections. 
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The joint boundary condition may, in general, be 
represented by a combination of an externa] moment, an angular spring 
and an imposed angular rotation, as shown at end A-of- Figure 5.4. A 
force boundary condition may be represented by ka = 0 and specifying 
a non zero value for My. A displacement boundary condition may be 
represented by setting Kn = © and specifying On. A spring boundary 
condition may be represented by setting On = 0, and specifying the 
appropriate values of kn and My: The boundary conditions at end B 
are similar, except that displacement in the axial direction is not 
restrained and the axial load P is applied at this joint. When 
necessary for numerical purposes, a linear horizontal spring is 
attached to joint B and the net axial force applied to the member is 


then P minus the spring force as explained in Section 4.2. 


A summary of the boundary conditions for the analyses 
to be presented subsequently is contained in Table 5.5. The symbol 
F, S or D has been appended to the model designation of Table 5.4 to 


indicate force, spring or displacement boundary conditions;,respectively. 


The material for all models is assumed to be elastic-perfectly plastic. 


5.3 Analysis of Member 3T for Joist AX05 

The center line geometry of this joist is shown in Figure 
5.3. The length of member 3T is seen to be 22.96 inches and the length 
of the joint member (between joints 4 and 5) is seen to be 2.36 inches. 
The axial force and joint displacements arising from a plane frame 
analysis of this joist, subjected to design loads, are shown in Table 5.1. 


The load displacement plots for the analysis which were run on this 
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member are shown in Figure 5.6. Each of these is described below. 


SiS «fl Force Boundary Condition Analysis for AX05 

The model of AX05 is detailed in Table 5.4. At design 
load the axial force in the member is 14730 lbs (Table 5.1). Table 
5.2 indicates that under design conditions the moments Mp and Mp 
acting on the member are -2.56 and -.116 in kips. The boundary 
conditions for the member are then those summarised in Table 5.5 for 
model AX05:F and a sketch of the loading condition is shown in Figure 
5.5(a). Applying these boundary forces in direct proportion to the axial 
load, a non-linear load deformation analysis may be carried out. The 
result is shown as curve AX05:F in Figure 5.6. The maximum carrying 
capacity of the member is 1.75 times the design load. The test value 
for this member, indicated in Table 5.3 and Figure 5.6, was 1.77. Thus, 
good correlation exists between the analysis with force boundary 


conditions, and the test value for this joist. 


bis. c Displacement Boundary Condition Analysis for AX05 


It can be argued that a force boundary condition on an 
isolated member will underestimate the collapse condition for the 
member as a part of a frame. The bases for this argument are: (a) the 
assumption that the moments applied to the member increase in proportion 
to the load Overestimates these moments if yielding occurs around 
the ends of the member, and (b) when the initial member becomes 


unstable, the end conditions tend to restrain deformations rather than 
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allow free rotation. To determine the influence of these effects, 
the member may be analysed by assuming the rotations at the ends of 
the member increase proportionately with load. This may be expected 
to produce an upper bound on the capacity of a member, whereas the 


force boundary condition may be considered to produce a lower bound. 


The member end displacements may be computed from the 
output of the elastic frame analysis of Table 5.1 and a sketch of 
deformations shown in Figure 5.5(b). The rotation of the member from 
the chord is obtained by substracting the chord rotation, obtained by 
dividing the relative vertical displacement of joints 6 and 5 by the 
member length, from the absolute rotations of Table 5.1. Detailed 
calculations are contained in Appendix E. The resulting rotations, at 
design load, are shown as the displacement boundary conditions for 
model AX05:D in Table 5.5. These rotations may now be increased in 
proportion to the axial load, and the resulting behavior of the 
member is shown as curve AX05:D in Figure 5.6. The member has a 


considerably higher failure load under these conditions. 


5.333 Spring Boundary Condition Analysis for AX05 


It can also be argued that the web members framing into 
a joint act as restraining angular springs on the rotation of the 
chord and a better model would be obtained by considering the top 
chord member to be elastically restrained at the ends. Analysis of this 
type was carried out by obtaining the stiffness of the angular spring 


in the following manner. 
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The elastic moments in members 4, 34 and 5 framing into 
joint 5 are shown in Table 5.2. If member 34 is considered to act as 
an angular spring restraining joint rotation, the spring stiffness 
may be obtained by dividing the moment of 2.205 in-kips in this 
member by the net rotation at the end of member 3T (6.058 radians 
as shown in Table 5.5), resulting in a spring stiffness of 367900 in- 
kips/radian. The external moment to be applied to this joint is then 
the moment arising in member 4, namely 4.789 in-kips. These values 
are shown as the spring boundary condition for AX05:S in Table 5.5. 
The conditions at end B are determined in a similar manner. A sketch 
of the model is shown in Figure 5.5(c). 

Increasing the forces in a proportional manner on the 
Spring boundary condition model of member 3T of AX05 results in the 
load deflection curve denoted as AX05:S in Figure 5.6. The behavior 


falls between that for force and displacement boundary conditions. 


oye et Analysis of AX05 Including Joint Length 


A final model of member 3T in joist AX05 was developed by 
including the length of the joint with the member. End A of the member 
was considered to be joint 4 of Figure 5.1, and was subjected to the 
force boundary condition arising from the sum of the moments in 


members 3 and 33. The boundary conditions are shown as those for 


AXU5J in Table 5.5. Joint 5 then became node 2 of the finite element 
model and the loads applied at this point, shown in Table 5.4, represent 


the summation of the components of forces in member 34 and one-half 
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the external panel point load. A sketch of the model is shown in 
Figure 5.5(d). The load deflection analysis for this model is 

shown as line AX05J:F in Figure 5.6. The critical load obtained from 
this model is at a load factor of 1.63 which is 7.9% less than the 


test value of 1.7/7. 


Some difficulty was encountered in this model in obtaining 
the unloading part of the load deflection curve. The curve is plotted 
in more detail in Figure 5.7. The reason for this behavior is shown 
in Figure 5.8 where it is apparent that once the critical load has 
been reached the deflections at the centre of the member begin to 
decrease because of large curvatures developed in the joint region. 
The load-deflection curve for horizontal displacement of end B is also 
Shown in Figure 5.6 which indicates an increasing displacement as 


the load decreases in the post-buckling region. 


Bot Review of Analysis for AX05 


The numerical values from which Figures 5.6 and 5.7 have 
been plotted are given in Tables 5.6, 5.7 and 5.8. The maximum load 
factors for each of the four analyses are underlined in Table 5.8. It 


is apparent that analysis AX05:F gives the best results for this member. 


5.4 Analysis of Member 3T for Joist AX0] 
The primary difference in loading conditions between 
joist AX05 and joist AXOl is that the external loads were applied 


directly to the top chord in joist AX01] rather than at the joints. The 
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joist geometry is detailed in Figure 5.2 and the elastic results of 
the plane frame stress analysis, for loads applied 2.915 and 4 inches 
from the centers of joints 5 and 6, respectively, are summarized in 


Tables 5.1 and 5.2. 


The loads acting on the model are shown in Table 5.4. 
Since the load adjacent to end B did not fall on a nodal point, the 
actual design load of 843 lbs was replaced by two statically equivalent 
nodal loads determined by the inverse lever rule as indicated in Table 
5.4. Boundary conditions were handled in the same way as for member 
AXO5 as described in Section 5.3. However, since the displacement 
boundary conditions give an unrealistically high collapse load, only 


force and spring boundary conditions have been considered. 


554 Force Boundary Condition Analysis of AXQ] 


The model of AX01 is detailed in Table 5.4. At design 
load the axial force in the member is 14660 lbs (Table 5.1). Table 


5.2 indicates that under design conditions the moments, My and M, acting 


B 
on the member are -.400 and -2.79 in-kips. The design load boundary 
conditions for the member are then those summarised in Table 5.5 for 
model AX01:F and a sketch of the loading condition is shown in Figure 
5.9(a). Multiplying these by a common load factor a non-linear load 
deformation analysis may be carried out. The result is shown as curve 


AXO1:F in Figure 5.10. Tris result predicts about 11% higher capacity 
than the test result indicated in Table 5.3 and Figure 5.6. 
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Sy gr Spring Boundary Condition Analysis of AX01 


The spring stiffness at joints 5 and 6 are obtained in 
a similar manner to those for AX05. The external moments to be 
applied to these joints are then the moments arising in members 4 and 
6, namely 2.604 and 2.973 in-kips. These values are shown as the 
Spring boundary conditions for AX01:S in Table 5.5. A sketch of 


the model is shown in Figure 5.9(b). 


Increasing the forces in a proportional manner on the 
Spring boundary condition model results in the load deflection curve 
denoted as AX01:S in Figure 5.10. As expected, this behavior predicts 


a higher capacity than the force boundary condition. 


54.3 Analysis of AXOl Including Joint Length 


Since both of the above models overestimate the strength 
of the member, a model including the length of joint with the member, 
similar to AX05J, was analyzed. The boundary conditions are shown as 
those for AXO1lJ in Table 5.5. Joint 5 became node 2 of the finite 
element model and the loads applied at this point, shown in Table 5.4, 
represent the summation of the components of forces in member 34 and one 
half the external panel point load. A sketch of the model is shown in 
Figure 5.9(c). The load deflection analysis for this model is shown 
as curve AXOlJ:F in Figure 5.10. The maximum carrying capacity of the 
member iS 1.76 times the design load. The test value for this member, 
indicated in Table 5.3 and Figure 5.10, was 1.64. Thus, this model over- 


estimates the capacity of the member to be 7.3%. 
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Di a Review of Analysis for AX0l 


The numerical values from which Figure 5.10 have been 
plotted are given in Table 5.9. The maximum load factors from the 
three analyses are underlined in Table 5.9. It is apparent that none 
of these analyses give very good results. All analyses predict a 
higher capacity than the test value. However, in contrast to AX05, 
the analysis including the joint yields a better value tnan that 


excluding the joint. 


5.5 Analysis of Member 3T for Joist AX02 


There is practically no difference between joist AXOl 
and AX02 except that joist AX02 was fabricated with larger joint 
eccentricities. External loads were directly applied on the top chord. 
The joist geometry is detailed in Figure 5.2 and the elastic results | 
of the plane frame stress analysis, for loads applied 2.515 and 4 inches 
from the centres of joints 5 and 6 on member 3T, are summarised in 


Tab ves 5 rand: x2’. 


The loads acting on the model are shown in Table 5.4. Since 
the interior loads did not fall on nodal points, the actual design load 
of 843 lbs was replaced by two statically equivalent nodal loads as 
indicated in Table 5.4. Boundary conditions were handled in the same 
way as described in Section 5.3. However, since displacement and 
spring boundary conditions overestimate the collapse load, only 
analyses for force boundary conditions, with and without the inclusion 


of the joint length, have been considered. 
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5 Decl Force Boundary Condition Analysis of AX02 


The model of AX02 is detailed in Table 5.4. At design 
load the axial force in the member is 14670 lb(Table 5.1). Table 5.2 
indicates that under design conditions the moments My and Mp acting 
on the member are -1.784 and -2.969 in-kips. The boundary conditions 
for the member are then those summarized in Table 5.5 as model AX02:F, 
and a sketch of the loading condition, at design load, is shown in 
Figure 5.11(a). Multiplying these by a common load factor a non- 
linear load deformation analyses may be carried out. The result is 
Shown as curve AX02:F in Figure 5.12. The maximum carrying capacity 
of member is 1.69 times the design load. The test value for this 
member, was 1.55. Thus, this analysis overestimates the member 


capacity. 


Sigs Analysis of AX02 Including Joint Length 


The joint model was developed similar to joist AX01. The 
boundary conditions are shown as those for AX02J in Table 5.5. The 
joint length is 2.97 inches and the length of each element is assumed 
as 1.5 inches. Hence the joint length is approximately two element 
lengths, so joint 5 becomes node 3 of finite element model. The loads 
applied at this point, shown in Table 5.4, represent the summation of 
the components of forces in member 34 and one half the external panel 
point load. A sketch of the model is shown in Figure 5.11 (b). The 
load deflection analysis for this model is shown as curve AX02J:F in 
Figure 5.12. The maximum carrying capacity of the member is 1.56 


times the design load. The test value for this member, indicated in 
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Table 5.3 and Figure 5.12, was 1.55. Thus, very good correlation 
exists between this analysis, with the effect of joint length, and 


the test value for this joist. 


ysis ae! Review of Analysis of AX02 


The numerical values from which Figure 5.12 has been 
plotted, and the load factors for these analyses, are given in Table 
5.10. It is apparent that analysis AX02J:F gives the best result for 


this member. 


5.6 Effect of Transverse Load Position and Yield Stress 


The influence of the location of the transverse loads 
on the member is illustrated for joist AXO1] in Figure 5.13. Two 
solutions were obtained: 

(a) With the transverse loads placed 6" from either end. 

(b) With the transverse loads placed 3" from either end. 
Force boundary conditions were used (AX01:F of Table 5.5). The 
numerical values for the analyses are shown in Table 5.11 and the 
ultimate load factors are 1.84 and 1.35 for cases (a) and (b), 


respectively. 


It is apparent that the positioning of the transverse 
loads on the span has a dramatic effect on the carrying capacity of the 


member. 


The influence of yield stress on the ultimate load factor 


is illustrated for joist AXOl in Figure 5.14. The numerical values 


ie tna’ -96e a = ie Rip 
x aes 
deeded bale a A a a) 


afdeT at nsvig 516 .2zo2ylens Szor} rot ue 
Not tfuesr ery. ond araticlly 4: — Ribicoetinn su a 


4) GORGES cmtae 


ebsol ‘oevayensad ond 20 ‘notssoet aig to OTE, ont sila 

owl .£1,2 gword at f0XA derot vot betsvdaulit 2f yednsm sad no” 

- sbentstda, gram sles 

.bns vorltis mort. "3 bsosigq ne ead oe ‘ae a eo 
Jone rortite not “€ ObdgTa) abso? ‘oovavaland lait AMM” ‘ 
aT . (2:2 efdsT Yo 4: fOXA) bead siaw ‘efiott tbntos Vabnuod eoted a2 

hs Sn 8 1c ie i i Si hin a 

+ (a) bri "Ca) 28289 dt BET bed hn 7 

eel o@! . leon deveafa OFT Be a oxcaee ate 


ne se i oy ohagt a pete. Pita tt 


10 REP aae 


a _ 
i 


a dan ve 
aii fone tm Zh amr 


eee ee eee a 


66 


for the analyses are tabulated in Table 5.]1. The analyses 

were carried out for loads positioned as in case (b) above, and result 
in ultimate load factors of 1.35 and 1.42 for FY = 60.3 Ksi and 

FY = 65 ksi, respectively. This result indicates that the ultimate 


load is also sensitive to the yield stress. 


5.7 Summary of Top Chord Analyses 
A summary of critical load factors for the different 
models that have been analyzed in this chapter is given in Table 5.12, 


where they are compared with the test results. 
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CHAPTER 6 
SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 
6.1 Summary 


The application of a set computer program for the large 
deformation analysis of inelastic beams has been investigated for a 
variety of problems. The programs were used to obtain maximum loading 
conditions for W8x31 sections which compare favorably with the inter- 
action curve of Galambos and Ketter. The inelastic program with 
initial stresses (PLAST]) has been modified ( to NLHAT) so that it 
is applicable to unsymmetric I sections with a residual stress pattern 
simulating that which could occur in a hat section. The hat section 
has been simulated as an unsymmetric I section and interaction curves 


have been generated for this type of section. 


The program NLHAT was then applied to models of the 
critical member in the top chord of three of the joists in Matiisen's 
test series in an attempt to determine if such a model can properly 
predict the failure conditions of the member when it is a part of a 
joist framework. The models considered the effect of boundary conditions 


and the effect of including the joint length as a part of the member. 


6.2 Conclusions 


The following conclusions may be drawn from the investigation. 
ee The finite element program used herein can properly account for 


the effect of large in-plane deformations. 
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ie When used with the augmented stiffness technique the programs 

can be applied to predict ultimate carrying capacities of members. 

Si The programs predict ultimate strength conditions for symmetric 
sections which are in good agreement with the beam-column interaction 
curves of Galambos and Ketter. 

4, The program NLHAT underestimates the fully plastic interaction 
conditions for unsymmetric sections and is, therefore, conservative. 

of The interaction curve for a beam-column with an unsymmetric 
section is considerably lower than that for a symmetric section when 
bending produces compression on the narrow flange. Thus, design rules 
formulated for symmetric sections may be inadequate when applied to 
unsymmetric sections. 

6. Models of top chord members which have either spring or 
displacement boundary conditions consistently overestimate the capacity 
of the member in a joist framework. 

ike The ratios of test to predicted load factors are 1.01, 0.89 and 
0.91,for force boundary conditions, for joists AX05, AX01 and Ax02, 
respectively. 

o. The ratios of test to predicted load factors are 1.09, 0.93 and 
0.99, for force boundary conditions applied to members augmented by 

the joint lengths, for joists AX05, AXO1 and AX02, respectively. 

9. The load factors for the members augmented with the joint lengths 
are consistently lower than those without, and this type of model gives 


the best correlation with the test values for the three joists examined. 
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10. Tne number of joists for which test values were available is 
insufficient to draw any firm conclusions about the reliability of 


the modeling techniques described herein. 


6.3 Recommendations and Observations 


The following recommendations and observations are made 
for future work. 
1k It is recommended that predictions of top chord capacities be 
carried out for other joists, associated with the Canadian Steel 
Industries Construction Council] Research Project No. 744, as test 
results become available. In this way the reliability of the 
predictive capacity of the program can be established. 
(be In view of the sensitivity of the predicted failure load to 
the location of transverse loads it is recommended that an experimental 
study of this effect be carried out. | 
8. In view of the sensitivity of the predicted failure load to 
yield stress it is recommended that improved techniques be developed 
for measuring the mechanical properties of hat-shaped sections. 
4, It is observed that the models of the members with force boundary 
conditions, but without joint lengths, considerably overestimated the 
capacity of two of the three joists available for study. This may be 
a result of the inadequacy of the boundary conditions imposed on the 
model. However, test observations indicate that there is considerable 
deformation of the cross-section prior to failure, which cannot be 


predicted with the present programs. It is possible, therefore, that 
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the discrepancy between predicted and observed collapse conditions 
results from the fact that distortion of the cross-section produces a 
decrease in carrying capacity of the section, and it is recommended 


that programs be developed which can adequately account Lor this-ettect. 
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Table 1.1 Section Properties 


Note: Units are in inches 
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Table 4.1] Comparison of Classical and Numerical Solutions 


for the "Elastica" 


DEFLECTION (inches) 


Paar | Wang's Result(22) 
| (Theoretical Solution) 
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Table 4.2 Inelastic Solution for a Beam-Column by Program PLAST 
{ 
Total Load | Net Load | Vertical Deflection *% Plastic Penetrat- 
in lbs. | ities in inch ion at x = 
200,000 | 200,000 | 0.2739 Nal 
747,200 | 253,115 0.3436 30 
747,500 | 253,104 0.3461 3] 
748,000 253,068 0.3506 33 
749,000 252,927 0.3611 | 36 
749,500 | 252,811 0.3672 37 
750,000 | 252,655 0.3742 39 
750,500 257,448 i 0.3822 | 4) 
| 751 ,000 252,184 0.3913 43 
| 785, 000 245,673 0.5447 | 62 
| 760,000 236,000 0.7434 ee i 
i 
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* 
~ Plastic Penetration = [d - yo 


- ¥3 Jx100/d (see Figure B-2) 
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Table 4.3 Moment-End Rotation for Beam Column” 
/M (1bs-in) | 6 
120,000 0.445 
128,400 0.508 
134,176 0.602 
129,703 0.795 
134,163 0.706 
134,867 0.665 
134,621 | 0.685 
134,836 | 0.669 
134,906 0.657 
134,912 | 0.653 


134,907 0.849 


“Example of Section 4.3: 2=138.8",~ = 40..— = 0.8>F = 33)ksi- 
r Py y 
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5~ Results from Galombos and Ketter (6) 


hi a ee RMR ts) cho. , lh My ml Deiat a 


0.75 | 0.65 | 0.56 | 0.47 0.38 | 0.29 
O067950-55 | 0.47 #| 0.39 | 0.32 | 0.23 
0.58 /0.46| 0.39 | 0.32 0.26 | 0.19 
0.49°100.87 1 0.31091) 0.25 | 0.20 | 0.15 
0.40 or | 0-23 | 0.19 0.15 0.11 
0.30 10720 | O.15-9)..0:13 10 | 0.08 
0.20 /0.11 0.08 | 0.07 | 0.05 | 0.04 
| | | 
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Table 4.5 Comparison of Compared Wo Values with Interaction 


Curve of Galambos and Ketter (6) 


0.2 


0.8 0.13 
COR ESN 


[] - without residual stress obtained from PLAST] 
() - Galambos and Ketter with residual stress 


Unbracketted Quantity - without residual stress obtatned from PLAST1 


Eat 


| ANS. 6) 


| 88.0. 


Table 4.6 
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Moment-Curvature (M- 


Ratio and Bending Compression in Wide Flange 
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Table 4.7 Moment-Curvature (M-¢) Relationship for Zero Slenderness 


Ratio and Bending Compression in Narrow Flange 
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Table 4.8 Finite Element Interaction Analysis for Simulated 
Hat Section. (Figures 4.6 and 4.7) 
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Table 4.9 Interaction Values for Simulated Hat Section 


(Wide Flange in Compression: epee: 
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Table 4.10 Interaction Values for Simulated Hat Section 


(Narrow Flange in Compression : Oe = 0) 
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Table 5.1 Design Load Results from Elastic Frame Analysis 
Joist Yield Strece,| fviatabaece 1 Joine Rotation ae 
(ksi) ; uf 3T(1bs) | (inches) 
5 -0.0110928 -0.4091547 
AX01 65 14660 
6 -0.001 7483 -0.5159207 
-0.0135107 -0.4730675 
AX02 65 14670 
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Table 5.2 Design Load Bending Moments from Elastic Frame An 
Member AXO] AXO5 
3 3.356 | 3.693 0.754 
| -3.065 -2.839 -0.755 
4 (Joint)! 4 8.801 | 7.658 
5 2.604 | 4.79 4.789 
-0.40 -1.784 -2.56 
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Table 5.3 Summary of Test Results 


Design Load (1bs) 1686 1686 1686 


Total “Utos} Load 2759 2769 2984 
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Table 5.5 Boundary Conditions for Top Chord Models 


Model > 3 


6x10 ky x10 opxl0> 
AX01:F Force -4.00 -2790 | 
AX01:S Spring +2604 339 .6 +2973 


Force — 


Spring 


Disp. 


AX05J:F | Force +7658 etree -116 - 0 


Note: Units are K-inches 
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Table 5.7 Load-Vertical Deflection (10 Values at Each Node 


for AX05J:F 


Load Factor | 


Node Noo | 1 1.63 1.53 1.60 

| 1 eae 0 0 0 

2 148 | .226 - 229 -.12 
3 271 506 49 267 
4 355 707 259 474 
5 4401 | .828 403 .607 
6 412 .872 481 .669 
7 .393 .844 5 . 666 
8 .346 754 466 .605 
9 .278 611 . 388 496 
10 194 429 wae 351 
11 ; +992 22 143 181 
2 | 0 0 0 481 

{ \ 


Note: All deflections are in inches 
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Table 5.8 Load-Deflection Relationships for AX05 


AX05:F AX05J:F AX05:S AX05:D 


Load Factor | u(107') U(1o 14 ptaad, |} UC07!) 
inches inches Factor} inches 
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Table 5.9 Load-Deflection Relationship for AX01 


_. AX0O1J0:F | 


Load =| u(io!) | Load uf i9.'),.1 
Factor inches Factor inches 


Saha: 
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Table 5.10 Load-Deflection Relationships for AX02 


AX020:F | AX02:F 


| u(1o-!) u(107!) 
inches inches 


Load Factor Load Factor 
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Table 5.11] Effect of Load Position and Yield Stress on Joist AX0l 


Deflection at mid span for AX01:F 


Load at 3" from joints! Load at 6" from Joints 


F = 60.3 

y Fy = 60.3 Fy = 65 
Load u(107!) Load u(107!) Load u(1o7!) 
Factor inches Factor inches Factor inches 
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Table 5.12 Summary of Critical Load Factors of Joists 


Joist Load Factor Load Factor Test Factor 
From Analysis From Test Predicted Factor 


AX05:F : Lez7 
AX05J:F 
AX05:S 
AX05:D 


AX01:S 
AX01:F 
AX01J0:F 


AX02:F 
AX02J:F 
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WEB TOP CHORD 


BOTTOM CHORD 


FIG.1-1 TYPICAL OPEN WEB STEEL JOIST 


FIG.L2. TYPE C - SHALLOW HAT SECTION 
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Fic.1.3. Equivalent Section. 
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Fig.2.1. Beam. Column Model with Fig.2.2. Beam- Column Model with 


End Moments. 


Transverse Loads. 
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Fig. 2.3. Bending of Rectangular Section. Fig.2.4. Inelastic Stress. 
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Fig. 2.6. Effective Modulus versus Total’ Strain diagram (2! 
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Fig.2.7, Moment versus Total Strain Curves for 


Rectangular _Section(21) 
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Element of Column Length 
(c) 


(d) Load Deflection Curve 


Fig.2.8. Details of Von-Karmans Analysis (15) 
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Fig.2.9. Stress - Strain Distributions for a Steel section. 
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Fig.2.10. Typical Progressive Yield Conditions for a 
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Fig.2.1l Sample Set of Assumed Yield Penetration conditions 


for Determination of Auxilliary Curves (11) 
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Fig. 2.12. Typical Set of Auxilliary Curves (11) 
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Fig. 2.14. Moment- Thrust Curvature Relationships for 8WE 3l (6 
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Fig.2.15 —TYPICAL NUMERICAL INTEGRATION PROCEDURE 
TO OBTAIN END SLOPE (6) 
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Fig.2.16 —TYPICAL MOMENT VERSUS END 
ROTATION CURVE (6) 
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Fig-2.17. Maximum Carrying Capacity Interaction Curves (6) 
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Fig.3:2 Basis Functions for Cubic Polynomial 
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Hand Computation 


Finite Element 


———— Symmetric Sections (6) 
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Fic. 4.7 Interaction Curves for HAT SECTION (Fy =33) 
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Fig. 5.5. Loading and Boundary Conditions for AXO5 Models. 
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Fic. 5.7. Load - Deflection Curve for AXO5vJ. 
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Fig. 5.9. Loading _and_ Boundary Conditions for AXO.1 
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Note: Units are inches and pounds. 


Fig.5.11. Loading and Boundary Conditions for AXO2 Models 
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APPENDIX A 
ASSUMED DISTRIBUTION OF RESIDUAL STRESSES 


The residual stress distribution for a typical wide 
flange section is as shown in Figure A.1. The value Of os in terms 


f-o is gi 
of oo.» 1S given as 


ubte,. 
ont bt + wid-2t) rc (A-1) 


where b and d are the width and depth of the section and t and w are 

the thickness of flange and web respectively. The distribution of 
residual stresses on an unsymmetric I-section equivalent to a hat shaped 
section (Type C, manufactured by Great West Steel) has been assumed as 
shown in Figure A.2. This figure indicates the residual stress in 

terms of three controlling values. (re. O_ and One): However, the 
two equations of equilibrium which these residual stresses must satisfy 


are constraint equations which reduce the number of independent values 


to one. In this Appendix the values of o 


bie and sino are evaluated as 


FUnCTIONS Of ose. 
rc 


Equilibrium of forces (that is, =M = 0 and =F=0), are 
considered, assuming the forces act at the middle of the fibre, and taking 
the moment about the middle fibre of wide flange. Considering the 
stress diagram and section properties of the unsymmetric section shown 
in Figure A-2 and Figure 1.3, and writting the equilibrium equations, 


one obtains; for the <F = 0; 
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(e) ano 
scat! es 2 i ee ae ak os 
1 "rc, ort "rey ort 
: 23> 3% rt. Ort 86 8. Sr ormee re) 0 (A-2) 
ae “rt °rc ort 2 


where (b, ,b,,b3) and (a) 285,43) are widths and depths of the narrow 

flange; web and wide flange respectively. on is the maximum compressive 
] 

stress at the top of the web and is considered to be uniform in the narrow 


flange. o., is the tensile stress at the bottom of the web and the middle 


Yat 
of the wide flange and a is the maximum compressive stress at the tip 
of the wide flange. Compression and tension are indicated as negative 
and positive, respectively. 


Simplyfying Equation A-2, it can be written as 


a re (244) ae (Ont pipet 2 (Chars = Ona 0 (A-3) 


Considering now the equilibrium of moments, the following equation is 


obtained. 
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Equation A-5 can also be written as 


2 
ata a,b.a Aad a,b,a a, b 
jag os AbD ag a Po. Dope te. 801 o 
pc, baby (apt mime mete cor gos tary saa Aime med 


Putting the value of oe in terms of Ong? from Equation (A-6), into 
1 


Equation (A-3) 


Substituting 
ab, aa a3 ay 
D = (abot aga) aaape Po” aj+ a3) + q+ a] - (2ajby+ andy )(q~ + =) 
yields 
a.b.o 6a.,b 
ga Ody 
ge an a Lappe io at a) + 3a4t 4a. (A-8) 


Putting the value of On, From Equation(A-8),into Equation (A-6), and 
simplifying, it becomes 
(3a 2a,.)a b.o 


2 AS eee atc 
"rc, = ay (A-9) 
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The values of O,, and On, are expressed in the program in terms of 


C 
] 
aoe by Equations (A-8) and (A-9). 
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APPENDIX B 
DETERMINATION OF CRITICAL STRAIN LOCATIONS 


The purpose of this appendix is to determine the points 
in the web of a simulated hat section, where changes in the modulus 


of elasticity occur. 


The equation of any linear variation of strain through the 


web may be written as 
= mix +e (B-1) 


where (x,e) are co-ordinates of any point on the line, m is slope of 
the line, and b is the intercept on the ec axis. Assuming linear 
distribution of residual strain, as shown in Figure B.1, and sub- 
stituting the co-ordinates of A and B into Equation B-1, the constants 
mandc can be obtained for the line A-B. Sign aeeiahione are also 


shown in Figure B-1 and Figure B.2. 
Evaluating Equation B-1 atA yields 
- re, = m(X, - ay) +b (B-2a) 


and evaluating Equation B-1 at B yields 


OF oes m[-(a, + a5) + XgJ +b (B-2b) 
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Substracting Equation B-2a from Equation B-2b gives 


“page mLxg - (a, + ao) - xg + ay] 


Hence 


pees mall (B-4) 


Substituting the value of m from Equation B-4 into Equation 


B-2, the value of b can be obtained as shown: 
(a, + ay - Xg) + b (B-5) 


By algebraic manipulation 


rc, 1X9 a (a,+ an) ] e Ent (Xg iG ay) 
b = Pi eae OE, Ee (B-6) 


The above values of m and b define the residual strains in the web. 
The critical strain points where changes in elasticity occur are shown 


in Figure B.2 and, if the strain from superimposed loads may be expressed 


as f + cx,are located from the following expressions. 
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at yy 

fir ic xt wet, by = ey, at Yo (B-7) 
~€y] at ¥3 
-€ aie Ya 


Solving Equation B-7 for each case, the strain locations are 


obtained as 


fest Do) 2) 5 
EN Loe Ls 
ae Coat (B-8) 
Wiel iy aces 
= Sse Se Bs x 
raid Catam (B-9) 
fib Wes 
i ee (B-10) 
ft+tbte 
P y2 
yee oo+ mt (B-11) 
Ween Xa ay (B-12) 


The strain locations y to Yq define the begining of positive and 
negative strain hardening from the centroid of the section. The 


distance from the centroid to the top of the web is denoted by Ys° 
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Fic.B.1. Residual Stress and Strain Distribution in the Web. 
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“*— Exact line Substitutes for 


Exact Quadratic Expression. 
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APPENDIX C 
HAND COMPUTATION OF INTERACTION CURVE (= = 0) 


In order to determine the fully plastic stress resultants 
for the given equivalent unsymmetrical I-section shown in Figure 1.3, 
it is first divided into small strips as shown in Figure C.1(a). The 
centre of gravity and the equal area axis of this section are at 0.453 
and .251 inches from the extreme fibre of the wide flange. The stress 
diagram for the fully plastified section is also illustrated in Figure 
C.1(b).The area of each of the strips is calculated and the moment of 
area of all these strips is evaluated about point O (Figure C.1(b)),and 
given in Table C-1. The location of the discontinuity in the yield 
stress (point A of Figure C.1(b) is assumed to move to the edge of each 
strip in succession and positive and negative plastic areas, the net 
area and the net moment of area are calculated as shown in Table C-2. 
The moment of area about the center of gravity is then obtained by 
adding the net moment of area about point 0 and the product of the area 
times 0.4536 inches (i.e. - the distance of the center of gravity from 
the bottom of wide flange). These computations are shown in Table C-3. 
The top of the fifth strip is the equal area axis and at this section, 
it is observed that the net area is equal to zero and the net moment of 
area therefore represents the fully plastic moment in pure bending. 
Areas at the other penetrations are divided by this net area to obtain 
py values and corresponding values of 77 are obtained by dividing the 
total moment of area at every section by this net moment of area. Detailed 


calculations to obtain these values are shown in Table C-3. 
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Table C-1 Calculation of Static Moments 
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Table C-2 Calculation of Net Static Moments 
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Table C-3 Calculation of Points for Interaction Curves 
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APPENDIX D 
PROGRAM NLHAT 


D-1 Preparation of Input Data for Program NLHAT 


Program Limitations 
Number of elements < 20 
Number of nodes < 20 
Input: 
(1) ALFA (1 card (Format F10.0) 
ALFA is a value used in numerical differentiation of the 
equilibrium equations (normally specified as -.05). 
(2) (NEL, MAX, NI, IPRINT) (1 card) (Format 418) 


NEL : Number of elements 


MAX : Maximum number of iterations in the Newton-Raphson 
procedure 
NI =: Number of subdivisions for the application of 


Simpson's rule (1<NI < 24) A good choice is NI=5, 
unless the elements are too long. 
IPRINT: Indicator for intermediate printing. The Jacobian 
matrix is specified at each step when IPRINT = 1 
(3) (XL, ALSPAZ§UA3,=B1, B25 BS) it licard) (Formatl7F0.0) 
XL ; Beam Length 
AlSMA2S (A3, BI (B23EB3s: cross sectional dimensions as 


shown in Figure 1.3. 
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SE We icewels, EY, CEY2 wSreymulecard) (Format 6F loud) 
ET], ET2, ET3 : Three moduli of elasticity involved 
in the trilinear constitutive equations. 
Eh, EY2 : The strains corresponding to changes in 
modulus of elasticity 
SRC -: The compressive initial residual stress 
whose distribution is shown in Figure A.2 
(5) INUMSP, EPS) (1 card) (Format 14, F16.0) 
NUMSP : Number of springs 
EPS : Specified accuracy for convergence criterion 
In case the number of springs is not zero then for each 
spring a separate card is included in sequential order. 
(6) (NPS(I), NDEG(I), ESP(I)) (one card for each spring)(Format 214, F16.0) 
NPS : Nodal point number at which a_ spring is attached 
NDEG : Local number of the degree of freedom (1,2,3, or 4 
for u, w, u' or w', respectively) 
ESP : Spring stiffness associated with the degree of 
freedom 


(7) IB(1), IB(2), IB(4*NEL+4)) (Format 1018) 


IB(1) =: O for an unconstrained degree of freedom 
1 for a constrained degree of freedom 


If IB(I1) = 1 then BND(I) is the value of the displacement for 


constrained (1th) 


(8) (BND(1), BND(2),...BND(4*NEL+4)) (Format 8F10.0) 


BND : The displacement specified for the degree of 


degree of freedom 


freedom 
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9). X(1)., X(2),..-X(NEL+T) (Format’ 8F10.0) 


h 


X(I) : X force (transverse) at the 1 node 


ieee 2 NELSL (Format or 1020) 


h 


LGU :, AXlal force at the i node 


11) XM(1), XM(2), XM(NEL+1) (Format 8F10.0) 


th 


XM(I) : Moment at the I~ node 


(12) NF, NF1 (Format 218) 
0 Initial trial vector on cards 
NFc= I initial tial vector in file 15 


? Initial trial vector in file 16 


Q last trial vector will be written on file 15, file 16 
NFl = will be left unchanged 
1 last trial vector will be written on file 15 but 
converged solution will be written on file 16, erasing 
previous vector 
(13) U(1),...U(4*NEL+4) (Format 8F10.0) 
These are the values of the initial trial vector and cards 
containing this are included in the data preparation if NF 
is specified as zero. 


D-2 Program Listing 


A listing of the program follows: 
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ao 


C MAIN PROGRAM FORK NCNLINGEAR ANALYSIS OF UNSYMMETRKIC 


did 


eo 


25 


bee? 


PMPLICIT REAL* 8 CA-H J0-2) 

DIMENSION U1100),D0U0100),4 1025) 
CGMMON/ POAD/Z(21) »XC21)—9XM(21) 
COMMUN/TBNO/I6 4100) 

COMMONS BNUCUN/ ENC( 100) 
COMMON/SONSTEGQ/ ETL ,ETZ, ETB, EYL, EYZ 
CUMMON/ SECT/SXG sAlL Ad 9A3 91982433 
CUMMUN/ RESTO/SRCULsSKT, SRC 
COUMMON/FINITE/E 

CUMMUN/NCHK/SIPRINT 

COMMON/ SPRING/ ESP140),NPS14C),NDEG(40) 
COUMMUN/DERIV/ALFA 

KEAD (5,100) ALFA 

WRITt(6,1234) ALFA 

READ{(5,20CU) NEL »MAKy,Ni yg LPRINT 
WRITE(6,iL1iad) NtbLyMAX NIT, LPRINT 
FORMAT (iHO,4110) 

READ(D,13C0) XLyAlegAcdsAS gtlytdy bs 
WRITE(O, 1 1e54) XLeAbpAdsASgbLdeBe BS 
FOURMAT(LHO,GE11~.5) 

ma Avtos LVCOO MET Ts ET2sETSs EYL EY!) . SRC 
Sete oc oo) VET ESET 2p ET Ss EYL SEY? SRC 


155 


I-GEAM 


D=((AL* Be tNs¥B 5) Fee FALHALTASIF(OCFALFEBIS(AL¥%B CC) tS oF ASH 4. FAL) —{ 
2eo®ALFBL FAL OC) ¥(5e*AStHtC 2 ¥AC)) 


SRO1=(A3*85*(35.¥*A34+22¥*AZ)/0)*SRC 


SRT=(AS*BS5¥*( (6 -*AL¥F OLS LAL 52) 2% (ee *ALHALHAS) 430 ¥AS44.*AZ)) *SRCE/D 


WRITE(O,1i234) SRC i, SRT 


AGH ae DALE( ACE EL) FU CSFLACHALIFA SEB AK (ACEU. DEI ABt4AL) ) I /S1( AL Bit AL¥B Zt 


=As*B 34 


WRITE(6,1235) XG 

FORMAT CULHO,1LUX g' XG=" ,Ei5.5) 
READ(D,12CU) NUMSP,EPS 

DQ 55> 1=1,49 

NPS{I)=0. 

NDEG(1)=0. 

ESP({1)=90. 


TF (NUMSP.2G6T.G) REAU(IS,icJi) (NPSTT) »NDEGTI),ESP(I),1=i,NUMSP) 
1200 FORMAT(14,F1E.C) 


FURMAT(214,F1l6.U) 
NPUTH=NELFi 
N=4*NPGOT 


H=XL/NEL 

READ(I5,2)00) (I5(1) ,l=iyN) 
WhRATG( oped lii) {lo(%),i=i,N) 
READ{>,1L 0000. (BNOU(1),1=1,N) 
WRITkE( 66,1254) {BNO(1I),1=1,N) 
KEAD(Dd,1900) (XC1T),1=LeyiNPOT) 


WRITE(6 p1le54) {X({1),1=1,NPOT) 
READ 55-1000) C2012 ,1=1,NPCGT) 
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WRITE(G,1234) (Z(1),1=1,NPCT) 
READL 5,2 9003 NF ,NFL 
IF(NF.GT.O) GO TO 234 
REAU (5, i900) {U(1),1=1,N) 
GO-TG 23-7 
254 [P)5=1 
NNK=L44+NF 
DOM ASD f=L,NecCT 
T1T1i=44#( 1-1) 
READ(NNK' I 15) (UCLTT+KR) sKR=L 94) 
£35 PRS=TP5+'1 
237 CONTINUE 
WRITE (6,123%) (U(1),1=i,N) 
CALL NEWRAPTIWN, LU, CU, NPOCT,cPS 9MAXgMI,NI) 
WRITE(6, 3600) {UC 1), f=1,N) 
30900 FURMATULHD LUA, RESULTS 4 y/S/ glLOA, TU LD X eg Wty LOK, tUPY gpLEXye* WP ',// y 
i (4620.7)) 
{FINFi.ebG.9) GOT Te 
ITlo=l 
Der 535 I=i,NPUOT 
Tli=4*(t-1) 
WRITEClot 116) (UCT IT+KR) »KR=1,4) 
335 TFO=1isti 
77 CONTINUE 
WRITiL(6,5000) 
5000 FURMATULHOg//// 4 6X, "ELEMENT NOS LOX, BENDING MCMENT* 5X, 
1 *MOu. AXTAL FORCE, BX," AXTAL FORLE*, 5X, 
Ps UY TEL POX MYT 24+! 5X, YT FS5—F* 9 5X, fT YITF4-+'/ 
3 TYW1t—',5X_* YW2t—* 9 5X_* YW3—+? gDXKy VY W4—-4+!* 5X, 
Remy Bree 5 OK ey YGF LA—F OKs YOR S—+8 5 OX, YY BG 4-4" ,//) 
VG 83 T=L,NEL 
CALL NMP{1,N1IT,U,N,1) 
683 CONTINUE 
TFINUMSPTEG.O}) GO TOIL? 
DU i190 I[=1,NUMS$? 
NN=NPS{T) 
N2=4*NN 
J=ENDEG(T) 
oJ is (90,597,955) _J 
96 EA=KANNIFESP 11) *UCN2-3) 
GU) Fe<99 
97 EX=ZINNI-ESPULD *UtN2-2) 
GB TO 99 
G8 EX=EXMOCNNI-ESP( LIFDATANZ(U(NZ—-1),01-4#U0N2))) 
99 WRITr(6,5500) NNoJgEX 
100 CONTINUE 
5500 FURMAT (LAD 77/7-,10X,*NEI FORCE AT NODE", 14,"DEGREE UF FREEVOM',14, 
i TE KIEL RNALEFORELY,ELZ.>) 
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DQ 20 T=1,NPCT 


19 
20 U(I4* 1-2) =U04*1-2)tH¥(I-1) 
WRITE(G,4000) (U04* 1-3) ,U(4*1-2),1=1,NPO0T) 
4000 FURMATILHC,///7/419K%_ COORDINATES S SS aL OXg SX LIX Z* pF J el et20e7) ) 
L000 FORMAT(8F10.0U) 
2000 FORMAT(1018) 
CAPE HEX TT 
END 
SUBROUTINE NMPCIT»NTyUysNyg IND) 
IMPLICIT REAL*8(A-H,D-Zj 
C IND=0 - NCO PRINT =i - PRINT UNLY =CTHER —- PRINT ALSO 


10 


C FOR 


DIMENSION U{N) 

DIMENSION WI{5C) 

DIMENSION YTF(5) ,¥Wi5) ,¥uF (5) 
COMMON/FAINITE/b 
COMMON/ANM/PN(4) »PPN(1454)9PM(4),PPM14,4) 
COMMON/SECT/XG , Ais Ad 9 AS 781 582583 
FiP(2)=6.%*Z*(Z2Z-i290)/H 
F2P(2)=3.¥*2¥*i-4.%*Zt+1. 
F3P(Z)=-F1P(2Z) 
F4P(Z)=3.2%*2¥*2-2 #2 
FiPP(2)=(12.*2-6.)/ (H¥H) 
F2PP(2)=(6.%/-4.)/H 
FBPP(Z)=-FIPP(2Z) 
F4PP(Z)=(6.¥%2-2e.)/H 

DO 10 IZ=1,4 

PN(1IZ)=0. 

PM(1IZ)=0. 

DG 10 JZ=1 94% 

PPN{I2,JZ)=0. 

PPM(I2Z,JZ)=0- 

CUNT INUE 

UL=U(4*1-3) 

UPL=U(4* 1-1) 

UR=U(4%1+1) 

UPR=U{4*1+3) 

WL=U(4*%1I-2) 

WPL=U14*I1) 

WR=U(4* 142) 

WPR=U(14*1 +4) 

SIMPSUN'S KULE 

NI2=2¥*NI+1 
TRINIT2sGEL38ANDNI2eLT.50) GO TO 13 
WRITE(6G,iiii) NI 
FORMAT(LHO,*ERKCK — NI=",I2) 
CREE SEX TFT 

WIC L)=Le/ Oe 

DM iD IW=L,NI 

W1I(c*IWw)=42/6.6 
wil(2*[wel)J=2./C. 
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2000 


300 


400 


WI(NIT2)=1.-/6.~ 

DO 4J0 J=iyNi2 

Z=DFLOAT(J-LI/SUFLUATINI2Z-1) 
UP=UL¥#¥FLIP{Z)4UPL¥F2P(Z)+UR¥*FS3P(ZIFUPR*F4P (2) 
UPP=UL¥*FIPP{(2) tUPL #F2PP( Z) 4UK ¥F SPP (Z2)4+UPK*F4PP(Z) 
WP=WLYFLP(ZItWELFF2P (2) tWK¥FS PU Z 2 tWPREF4OP( 2) 
WPP=WL*FIPP(Z)4WPL¥F2PPIZ) tWR*FSPP(ZIFWPR¥F4PPE2L) 
E=WP+0 .5*(WP*WEtUP*UP ) 

C=UPP*(1.+WwP)-UP*WPP 
C=C/((2.* E41.) 20SQKT(2.%E+1.).) 

CALL FORMOM(E, CaXNeAMeYTF a YWo YOSEF) 

[Ir (IND.EQ.0) 6C TG 300 

AXM=XM*(2e*E+L 2) 

XAN=VSURT(2.* C41.) FKNG2 OF ( 2 KEL.) FC EXM 
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WRITE(6,2000) Tg JSyXXM_XNy XXNy CYTE (LYTF) » IYTE= 154)% 


LOUYW(1Yw) sTYwW=i 4) s(YERF(IYBE), TYBFH1L 14) 
FURMAT(LH 9215 ,1L OX, SELB. 7p 4EL1 2 5/3EL1 2 5) 

If (INDLEQ.1) CO TU 400 
PN(LI=PNCL)I WIC J2*XN*FLP(Z)/NI 

PN( 2) =PN( 2) tWI (J) *AN¥FZP(Z)/NI 
PN(4)=PN(4) WIC 2 *XN¥F4P(Z)/N1 
PM(L)J=PM(1L) +WwI( J)*AM*FLPPIZ)/NI 
PM(2)=PM(2) 4+WI1 (J) *XM*F2PP(Z)/NI 

PM(4)=PM(4) 4WI(S)*XM¥F4GPP(Z)/NI 
PPN(L,LI=PPNCL pLIFWI(S)FXNFFLIP (ZIFF LPC LZI/NI 
PPN{1,2)=PPNI 1,2) WI (JS) *XN¥FLP(Z)*F2P0Z)/NI 
PPN(254d=PPN(i 14) tw1(d)FXNFFLP UL) SF4P(Z)/N1 
PPN( 292) =PPN( co ZdtWI(SVEXNFFZPIUZIRFSZP(ZI/NI 
PPN( 2,4) =PPN(254) 4+W1ITI) *¥XN¥F2P(L) *F4P(Z)/Ni 
PPN (4,4) =PPN(444) 4WI(5 2D¥AKAN¥F4PUZ)*FGPIZ)/NI 
PPM( ig id =PPM( ip LItWI(S)*AM*®FLPP(Z) ¥FLP(Z)/ NI 
PPMULy2Z)=PPM(iyg dd tW1( JS) FXMFFLPP(Z)¥*F2P(Z)/NI 
PPM( 1,4) =PPMi1 »,4) WI (J2*AMFF LEP Z)¥*F4P(2)/NI 
PPM( 2,1) =PPM( 2,1) wl (JS )*XM*FFOPP(Z)*FLP(Z)/NI 
PPM( dy 2) =PPM( 2,2 )4tWIT (SD FAMRF EPP (2) *F 2P(Z)/NI 
PPM(2 94) =PPM( 294) WILD E*¥XM¥FEPP(L)*F4GP(Z)/NI 
PPM( 4,1) =PPM( 4,1) +WwIl(J)*XM*FGPP(Z)*FLP(Z)/NI 
PPM(4,2)=PPM(4,2)4W1I(S )*XM*FFGPP(Z)*F2P(Z)/NI 
PPM( 494) =PPM(4 ,4) W115 2 *AMSFSEPP(Z)*F4P(2Z)/NI 
CUNT INUE 

IF( INCw~EQ.1) RETUKN 

PN(3)=-PNO1) 

PM(3)=-PM{(1) 

PPN( 1,3) =-PPN(U1,1) 

PPN( 2,5) =-PPN( isd) 

PPN(393)=PPN(01,1) 

PPN( 3,4) =-PPK(1,4) 

PPM 1,2) =-PPM( ls Ld 

PPM(243)=-PPM(2,1) 

PpM4(3,1)=-PPM(1,1) 
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PPM( 3,2) =-PPM{ 1,2) 
PPM(395)=PPM{1i,1) 
PPM(394)=—-PPM(1,4) 

PPN(2cei)V =PPN( i, 2) 

PP (5e9i)=PPN{1,5) 
PPN(3,2)=PPN12,3) 
PPN(4,1)=PPN(1,4) 
PPN{4_,2)=PPN(254) 
PPN(4,53)=PPN(13,4) 

RETURN 

END 

SUBRUUTING FORMUM( Es, Ce SNy SM ys YTF YW, YBF) 
IMPLICIT REAL*S (1 A-H,0-Z) 
COMMUN/RESITD/SKC1,SKT,SRC 
CUMMON/SECT/XG, Als Ac sASs6B1,B2,33 
COMMUN/CNSTEG/ ETA, ,ETZ es ETS sc YL, EYZ 
DIMe NSTON YTF (5) »¥wl 5) ,YoF15),STR{I 5D) 
SN=C. 

SM=0. 

DTD {=1,5 

YTF(Ij)=0, 

Yw{1)=0. 

YoF(1I)=0. 

ERC=SRC/ET1 

SRCUSSREA/ETI 

tRT=SRT/ET1 

STRAZI=HETI*EY1L 
STRIL)=STRIc)+ET2*(EY2-EY1) 
STR(3)=-STR( 2) 

STR(4)=-STRI1) 

GO TG 200 


THIS 1S A SGLUTION OF HAT SHAPED SECTION WITH CONSTANT C RKESI 


STR AT TOP 
WRITE(6,1100) 


FORMAT (i1HC,10X% ¢* ERROR -RESTOUAL STRESSES MUST BE INPUT POSITIVE") 


CALE AEXAT 

FLANGE 

SLOPEL=2.*1(ERTI+ERC)/B1 
X=XG-Al/2é. 

EX=E—?¢ eFC EX¥( Le ELEX ID F( E40. 5) 
YTF(L)=1 ES tERTL-EY2)/SLUPE1 
YTF (2) =(EGtckTI-EYL)I/SLUPEL 
YTF(S) = EBERT AtEYL)I/SLOPEL 
YTF (4) =(EBtERTL#EY2)9/SLUPE4 
Y¥iR(5)=81/2. 

STR1{5) =STRIGI + ETZ*(EB-ERCTHEY2) 
SU=STR(L)FETS* (EBHERTI-tY2) 
SUM=G-6 
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HO =S.04+8=] 75 

TFUYTF(1)) 45,45,25 
IF{YTF(T)-(81/2.)) 40,30,3C 
SB=SU+(STR{T)-SU)*(B81/2.-AC)/(YTF(1)-AG) 
SUM=SUM4#90.5*(564568) *({B1/2.-A0) 
GO=<EERLOG 

SUM=SUM#0.5*(SC+STR{1)) *1YTE(I)-AG) 
ADHVY.TF( 1) 

SO=STRC1) 

Gii-12 +50 
SG=STK(II-(USTR (14i2)-STR( DT) FYTF CLI /S(YTFCI4i)-YTF(1)) 
CONTINUE . 
SN=SN4+20*SUM*AL 

SM=SMA-Z. *SUM*A LEX 

GH(E6+259 

X=KG-Ali/2. 

EG=E—e oe FCFKF( 1 KOS *CHEX)D FC E40. S)-ERCL 
ERCCBeEEYZ) 2ils2idsZ2le 
STF=STR(4)4+ET3*(EBt+EY2) 

GO TC 239 

IFieeteEyi) 214,214,216 
STF=S5TR (3) 4+ET2*( EG+EY1) 

GO TC 239 

SFREGB-EY UT 2h e seth, 2Z29 

StFHSTEX EB 

GO TO 239 

bF(CEB-EYZ) L2La9e24 9224 
STFHSTRKRI 2) +c7T24(kEB-EYL) 

GO TG 239 

STF=STR(1LIFETS* 1 EB-EY2) 
SN=SNtUL*¥AL*ST SF 

SM=SM-BL*XAL*®STF RX 

CUNT INUE 


C 3OTTCM FLANGE 


2i 


20 


EPL eERI 1U,2Uil,2i 
SLOPEZ=2.*(ERT+ERC)/B3 
X=—-(ALF+A2+A3/ 2 --XG) 
EB=E—2.*C¥EX¥*¥(1.2.-0.5*¥C¥X)*¥(E40.5) 
YBFCULJ=HLlEBT+ERT-EY2)/SLOPE2 

YBF (2)=CEBt+ERT-~EYL) /SLOPE2 
YBF(Z)=(EBtERTHE YL I/SSLUPEZ2 
YBF(4)=(EBTERTHEY2)/SLGPE2 
YBF(5)=B3/2. 
STR(5)=STRI4) + ETS*(EBKERU+CY2) 
SU=STARILIFETS FC EBFERT-EY 2) 
SUM=Ve 

AU=0 2 

Di ( > 4: [=1,5 

lf (yer 1)) 46 546,26 

LEON BE LLFR LB S420 4 41,311,531 
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EBV LONI +1 

(Oxn0 Nene 

[2 .0+5 ek UERDeA, One L)# 
$2902 A S¥S—T sae ; 
S390 I2\ 01 ¥9-TAS4 sine 
$8932 \4 i¥ 3+ T wae 
SAMDJeNC S36 TaaeBas 
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(SYSeuaaquguee ananeel 2 =e Ws 
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aa 


310 
305 


325 


9B=504(STR¢C1)-—30)4(18687/2.-AU )71YB8F(1)-A0) 
SUM=SUM+0.5¥*{(SU4+SB)¥*(Bs/2.-AU) 


GOr TOI 101 


SUM=SUMtCO.5¥*(SC4STRI1))*(Y 6F(1)-AG) 


AG=YbFIT) 
SO=STR1(1) 
Ge Tt On ba 
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SU=STRI1)—-(STRCI4+L)-STRUI) O*¥YBF CE) /(YBF(I+1)-YBFE(1)) 


CONTINUE 
SN=SNtc 2 *¥SUMHAS 
SM=SM— 20 ¥SUM¥FA 5% X 
OONT DS 251 

K=— (AL+AL4+A3/2 -XG) 


PPA EBTEY 27 2 24%, 2415213 
STF=HSTR(I4)+ ETS + (EBTEY2) 
GO TO 258 

LPC EBPEY 19% ed, 2157217 
STFHSTRAISI+ET2*(EBTEYL) 
S04 TOeZBs 

ISEB EVENT O21 9 52195221 
SPS TI*Es 

GO TU 238 

PrP ESEY 2) * 228522395225 
SGrS=OTRAZ)+ET24(EB-EY1) 
ob TO 238 
STFHSTR(LIFETS F(t68-EY2) 
SN=SNt65¥* AS*STE 
SM=SM—-B3B*ASHSTE *¥X 
CONTINUE 


SIGN=1.0. | 
A=+(ERTHERCL)/ 42 


S= (CEKCL¥*¥(XG-(A1LFA2Z) PFERTEIXC-AL)JISAZ 


[F{CtA) 300,500,305 
SIGN=-1.0 

DOMS10/ UK= 14 
STR(K)=-STRIK) 


Yw(l)d=(EtB-SIGN*¥EY2Z)/(C+A) 
YW 2) =(EtB-SIGN*¥EYL) /(C+A) 
YW(3)=(E+Bt+STGN¥EVYL)/ (CHA) 
YW(4)=(EtB+SIGN¥EY?Z)/(C+A) 


Yw(5)=(XG-Al) 


STRISI=HASTRO4I + ETS 470E4+B-(C +A) *YW(5) +S TGN*EY 2) 
SO=STRK( i) teTS* (E+bt(C +A) *¥YW(5)-STIGN¥EY 2) 


AU=XG-(AL+FHA2) 
DOU 550 I=1,5 


DEM VWOL )-AOR s 03455345,325 
IF(YWCDI-YW(05)) 340,350,350 
330 SB=SCe(STRIT)-SO)*(¥W15)-AC)/ (y¥W1T)-AQ) 
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((1)48Y- ~(4 91) 48¥)\410 98" 


Bp ens es 
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Uvd-8a 65) faa 
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debit nny pelle 
(asa) At Lyasaatenaae 
(AFDINGSYSSADDE +8: 
AL a~a% 
(S ¥38N01 e+ (e LYS (AED SEHD CET OGTR 
he cence eee iTe=f 
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220 
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SN=SN40.54(SG456) *#(¥W(5)-ACI*B2 
SM=SM-(YW(5)-AO ) ¥52*150%*0. 54(YW(5)4A0)+156-S0)*0.5%*(A0t2.* 


—(¥W(5)-AD)/3.)) 


ot 


GU TO 400 
SN=SN#tC 25F(S5C#5TRC1) 0% (YW 1)-AvU ) ¥82 
SM=SM-(YWUL)-AC) Bc * (S0*0.5*(YW(1) 4AD) 41ST R(T )-SO) #0. 5*1{ AGF 


P2e*(YW(1)-AG)/3.)) 


AG=Yw(I) 

SO=STR{(1) 

SOLTONSS50 
SO=STR(I)-(STRIT4#LI-STRIT) 2 F (YW 19-AD)/CYW(T41)-YWil)) 
CUNTINUE 

RETURN 

PPteACBFEY 2d 51104 520% 512 
STFHSTR(G4) tt Ta*(EtB+EY2) 

GO T6539 

HPCE +E Y 1) 61 5145514, 516 

Se royed SPECT (et +EY i) 

GO 70 646 

LOLEAB ea) 8 O18 51S 49520 

STFH=ETL¥(E+d) 

60 TG 539 

IFN es —EYV29 B22 39922 ¢D24 

SHB SST RAR +b 2*1E+tG-E YL) 

GU 10 539 

SHfESTRAID + EISe E+ B=EV 2) 

SN=SN4+STP eB 2¥ AS 

RETURN 

cND 

IMPLICIT REAL*81A-H,O-Z) 

DIMENSIUN UN) 

COMMUN/FINITE/H 
COUMMUN/XNM/PN14) ePPN1454)_5 PM14),PPM(4,4) 
COMMUN/SOXNM/DPRH (495) pUPPN(4 9498) ,DPM(453),UPPM( 494,50) 
COMMON/CNSTEQ/ ETL, £LT2Z,E1TS,EYL,EY2 
COMMUN/DERIV/ ALFA 

DIMENSIGN APN14) yAPPN14 94) > APM(4) 5, AFPPM(4,4) 
DIMENSION BPN(4) ¢BPPN(4,4),8PM(4),8PPM(4,4) 
EY=EY1 

DU 50 K=194 

APN{K)=PN(K) 

APM(KI=PM(K) 

DU 50 L=l 4 

APPN{K,L)=PPNUK,L) 

APPM(K,L)=PPM(K,L) 

CONTINUE 

i) 100 J=1,3 

DELTA=ALFA*EY ¥*F 

LF (0 J—0 05-1914) *4)..6GT ed) De LTA=ALFA*EY 
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JJ=4*(1-1L) 4+) 
A=U(JJ) 
C ALFAeLT.0 ALTS AS INDICATOR TO USe A FIRST-GRDER DERIVATIVE 
LEGALFA) 5:005.:60,0 5700 
600 WRITE{(60,1000) 
1000 FURMAT(1LHO,*ERRGR. ALFA=0.") 
CALL. EX LT 
900 <YEbeRDEL TA 
U( JS =A-DEL 
CALL NMP (I,NIl alent) 
UL JJ) =A 
UG 559 K=194 
DPNUKsJ)I={APNCK)-PNIK))/JDEL 
OPM( Ky JJ=(APMIK J -—PMIKIQIV/JUDEL 
DU 550 L=1,4 
DPPN(KsL eg JJ=LAFPANIK yLI-—PPNUK,LDI/0DEL 
DPPMK,yL,JI=LTAEPM(K,L)-FPM (KL) /JDEL 
550 CUNTINUE 
GO af Os dao 
TOG U{JJ) =A-GELTA 
CALL NMPCIT,N1,UgN,O) 
DO 150 I1C=i,4 
BPN(IC}J=PN{(IC) 
bBPM{1TC)=PM( IC) 
DO 450 JC=1 5% 
BPPNTIC,JCI=PPN(IC,JC) 
BPPM{iC,JC)=PPM(1C,JC) 
150 CONTINUE 
U( JJ) =AtDELTA 
CALL NMP C1lyNIyl,N,9) 
Ui JJI=A 
DU 2uVd K=1L_4 
OPN(KyJI=LPNIK I-BPN(K))/(02-40DtELTA) 
DPM K,J)=(PM(K2-BPM(K))/(2.2*DELTA) 
DU 2u0 L=1,4 
DPPN(K,L es JJ =(PENIK,LI-BPPNIK,LII/(2.*DELTA) 
UPPM(Kyl,JI=(PEFM(K,LI-BPPMIK,L)I/S(2.¥*DELTA) 
200° CONTINUE 
100 CUNTINUE 
DG; 300 K=L 9% 
PN(K)=APN(K) 
PM{K)=APM(K) 
DG 300 L=ls4 
PPN{K,LJ=APPN{K,L) 
PPM(K,L)=APPM(KyL) 
300 CONTINUE 
Ke TURN 
END 
SUBROUTINE ABCGIN,U,NI) 
IMPLICIT KEAL*E(A-H,U-Z) 
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DIMENSION UIN) 

DIMENSTUN (494) 9B14 94) A494) 5614) 
DIMENSION APN(4),AP4(4),APPN(494),APPM(4,4) 
DIMENS TON ADPN (496),ADPM (4,3), ADPPN(4_ 493) y ADPPM(4, 4,35) 
CUMMUNSKXNM/SPNI4) pPPN 14 94) 9 PM14),PPM(4,4) 
CUMMUN/S DXNM/DPN(4,5 3) yDPPN( 4945 8)40D0PM( 458) ,.DPPM(4 54,8) 
COMMON/FINITE/S+ 

COMMUN/NCHK/IPRINT 

COMMONS POAD/S ZL eidy Xl 2b), xX4 (21) 
COMMON/TEND/16 (100) 
COMMUN/ SE NGUON/BND(1G0) 

COMMOGN/ SPRING/LESP140 ),NPS140),NDEG(40) 
Tii=i 

Tlé=1 

NP=N/4 

KOGUNT 2=i 

KUUNT=1 

DU 500 T=1,NP 
FI=VDATANZ(U(4* 1-4) ,(12+U(4%41))) 
CR=OC0S:2F EF) 

SF=US INCI FI) 

G(1)=-X 1) 

Ole j=-2 (1) 

G(3j9=-XM(1)*(1.24U(4*I)) 

(4) =AM(1)*U(4*I-1) 

PPE CEQ SNP) CU TCG 109 

CALL NMP(I,NI,U,N,02 

CALL DNMP (i,yU,N_NT) 

G6(i1)=G( i) tH*¥PM (4) 

Gl2)=G(2)+H*¥PNC1) 

G(3)=b513)4H*PM (2) 

G(4)=G6(4) FH ¥PN( 2) 

IPL FSEG th) GO TO itv 

GU LJ=G1 LP FH*APM( 3) 

G(2)=G6l2 PFHFAPHNI3) 

G{5)=G15) +H*¥APM(4) 

614)4=6(4) 4H+4PN( 4) 

CONTINUE 

D0 290 K=lyz 

DG 200 L=194 

L1=4*14+2*L 

TF(T.EQ.NP) CD TB SVZo0 


GA Z¥K-LI=EC( 24K -1)tHF(PPN(K,LI*FLCTI-5)4+(0FPM(K,L)-PPAM(L,K)) *U(TI-4) ) 
G(2*K)=G(2Z¥*K) +F ¥(PPN(K,LI*¥U(TI1-4)+(PPM(L,K J-—PPM(K,yL) *¥U(ITI-5)) 


IFtiseg@si) GOSTORZOO 


GU Z*EK-L)=C(2*K— 4 FHL APPN( K42_Ld*U(TI-9) + C APPM{ KtegL)-APPMIL,K+2)) 


i *y(11-8)) 


GI 2*K)=G6(0 24K) +h ¥(APPNIK+2,L)* U1 TI-8)4+0APPM(L,K+2)-APPM(Kt2,L)) * 


i UL II-9)) 


200 CUNTINUE 
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190 LFCI.NE-NPS{KULNT)) GO TO 2C2 
LL = NDEG{KCUNT) 
GO #titlb2Zlg2z2ida7rvF itt 
221 G(1)=601) + SP(KGLNT) *U14*1-3) 
GU TC 224 
222 Gla) =G( 2) +E SPI KCUNT) ¥U(4¥%1-2) 
GO T0 224 
22d GIS) =GI(B)FESPLKGUNT )¥F1¥*(1.24U(441)) 
614) =6(14)-ESP{(KGUNT )¥*FI*U{(44I-1) 
224 KCUNT=KOUNT+1 
60 70 190 
202 CONTINUE 
DY 210 K=1,4 
L=4-K 
GIK)=G(K)401.-1514*% I-L) )4+17U14*%1-LI-—6ND(4*I1-L)) ¥16(44%1-L) 
CUNT INUE 
vG 300 K=1,2 
DU 200 L=l,4 
IND=L/2-(L—-Llj/c 
LL=({L+i}/2 
ifdi pho oi GO, TOTS 10 
C1 2*K-1,L)=H*¥( ADPM{K42,L I+ APPN(K42,LL)*( 1-INDJ+t1(APPA(K4+Z,LL)- 
; KPPM(LL,K4+2))*IND) 
GCI 2Z*¥KyLJ=FF(ADPNIK4Z, LI TAPPNIK4tc,LLJ*INUF{ APPAILL sK4+c) -APPMIKt2, 
i 1L))*({1-1ND)) 
310 B642%K-i,L)=C. 
BI c*KyLJ=0. 
rid sti, Wr 9 Get TB! 320 
534 2¥K-1L,yL)=HF( OPM KR yLIFPPN{(K,yLLI*¥11L-INDI+1PPM(K,LLI-PPM{LL,K)) * 
L IND) 
B{(2* Ky LI H=Ee (UPI K,yLIFPPNiK, LL) FI NUt (PPMILL »K)-PPM{(K,LL) )¥(i-IND) ) 
(Hse st) Su Fo 4330 
429 BIZ¥K—-4L LJ Hb12*K—-igLk )FH*¥LADPM(Kt2ght+4) 4+ APPN(K42Z gL Lt+2)*(1L-INDI+ 
1 {APRM(K4ZyLLt+2)-APPAM(LL+2,K+2))*IND) 
B425KyL)=B( 24K gL) 4+H* (ADPNIK42,L44) 4+APPN(Kt2,LL4+2)*INOt 
lL {APPM{LL4+2,K+4+2)-APPM(K4+2,LL4+2))*11-IND) ) 
330 IF(1-EU.NP) GC TG 3440 
A( 2*¥KR-1,L )=H*¥(LPM(K,L4+4)4PPNIKyLL4Ec)¥(1L-IND) + (PPM(K,LL+Ec)- 


fo 
~ 
ee) 


1 PPM(LL4+2,4))*¥IND) 
AG 2¥* Ky LJ =H (OPAL KL 4+4)4+PPRNI Ky LLt¥c2) *INDt (PPNILL+Z2,K)- 
1 PPM(K,yLL4+2))*(1-IND)) 


340 CONTINUE 
DG 300 M=1;%4 
11=4%1+2*™ 
fF (L-EGer) GG TG 35d 
C(2*K-Llyl =U 24K-1,L) 40% (ADPPN(KF+Z)MyLI*UC T1-9) + 
i {ADPPM(Ktc yMyl )-AUPPM(M,K42,1))¥*U( TI-8)) 
C(2Z*K gL) =C(2*#K gL) +H¥(ADPPN(K+2 9M,L)*U 11-3) + 
i (AUPPM(MyKt2,L)-ADPPMI{K4+2)MyL))*U( 11-9) ) 
350 I[FLleEWeNP) GL TO 460 
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BI 2*K-id,yL=5(2*K—-LygLIFH*(OPPNU Ky MyLI*FUCTI-5)+ 


i (DPPM{(KyM,L)I-DPPM(M,K,L))*ULTI-4)) 
Ble*K,L)=B12*K ,L J tH* (DPPNIKM,L)*UCTI-4) + 
x (DPPM(MsKyLJ—-DOPPMIK yg My,L)I*ULTI-5)) 


IF(T.EQe1) GU TO 37o 
360 B(2*K—-1 pL) HB (2 FK—-1L yLI+He (ALPPNIK +2, MeL +4) %U( 11-9) + 
i (ADPPM{K 429M _L44) -ADPPM (BR ,K+t2,L+4) y*xU(TI-8)) 
b( 2% KylLI=B (eK aL DFH®(ADPPNUK4E 2 yM,Le44) FU( TIT -&)+ 
1 {ADPPM(M,K 42 9144 )-AUPPM (K+e9M,L4+4))4%0( T1-9)) 
370 CONTYNUE 
AC Z*®K-isl =Ale *¥K-1_,L) tH*¥IDPPN(KM,L4+4) *U(TI-5)+ 
ip (OPPM(KyM,L+4)-DPPM(M,K,L+4))*U(11-4)) 
AL2™H gL J=A(2¥*K gL PHH*(DPPN(K M,L44) FULL 1-4) + 
A (DPPM(M,KyL +4) -UPPM(K »M,L4+4))*U0(TI-5)) 
aVU0 CONTINUE 
D(oa_94)=b03,4)-XM(70) 
b(495)=3(14,3)4dAM (1) 
DO 450 L=l1y,4 
[L=4*{[+L-4 
OO 400 M=1 94 
[FCT Eg. NP] GG TOG 410u 
A(LsMI=ALoM)¥(I-LeUCIL)) 
419 6(i,MJ=B(LM)¥*¥CL-IBLIL)) 
Peer see.) B(i,MJ=o(iy»M)+1801L) 
PRCT .EQe1) Gt TO 400 
Clty MPHC OE, M8 OPT Bit Me) 
400 CUNTINUE 
422 ITFCLTNESNPSUKOUNT2Z J) GU Tu (421 
KK = NDEG(KUUNT2Z) 
GU Ti) (42974257424), KK 
429 B(1,id=B0L,L)+eESP(KLUUNT2) 
GO TG *425 
425 8l292)=8b (e292) +E SP(KOUNTE) 
GU TG 429 
424 63, 3)=0(03.3)+CSP(KGUNT 2) *CE*CE 
6 (494)=614,4)4+LSP(KUUNT 2c )¥SFR*SF 
843.4) =0(53.4) +E SPUKUCUNT2Z )*(-SFeCE+E 1) 
3(453)=B (74, 3)+ ESP (KOUNT 2)* (—-SF*CF-F I) 
429 KOUNTZ=KOUNT?EH] 
60 Tu 422 
421 CUNTINUE 
Trtid wos i? GC TO 450 
WRITE(aae' Tid) ((C{1T Py JP), JSPHi pt) yy TPHlL,4) 
14224 i2t1 
IFLIPRINT.EGe1l) = WRITEL(6,1313) ((C (IP, ISP), IP=ie4),1P=i94) 
WwRITTeE(1l2'tic) ({5(1P.JP) »JP=1,4),1P=1 54%) 
1Aceqi¢+. 
IF (IPRINTAEWe1L) WKITec(o, 1313) ((6(IPyJP),IP=in 4), [1P=194) 
LFLASEUL Ne Gu TU 460 
wkTTée(le* Tle) ((ACT Pe JP) eo IPHl ye 4d ep TPHly 4) 


fs 
Oi 
CG 


ee 


em ot | 
wa (as ¢ ~ anh 
gag 250 ee ae ee 
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ee 
* ig : poor 
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CL 1+3 412-18 (anette bbe 
+ et LaAete= 1=(STHUOAIS2 3 +66 50 oh ube | 
; lene (33 A= =sth Ww 


| 98 OT a Erae) 
(#ei@digtes i= =F ¢ fhe VL) hd) eee 
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469 
461 


losis 


490 


460 
500 
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iS 


30 


1 


lig=Jlig+l 

PEA PRINGT SEQ. 1) 
DU +61 IP=i9% 
G{TP)=-G(1P) 
WRAFtdi ld? lid) 
file=iids+ 
WRITrt(o,1513) 
FURMAT(iK »4E20.7) 

DU 460 K=1 44 
APN(KJ=PNIK) 
APM{KJ=PM(K) 

DU 490 L=1,%4 

APPN(K yi J =PPN( K re Pe 

APP M(K,LJ=PPM(K,L) 

OG 490 M=1,u 

ADPPN( Ky Le MI=DEPN( Kyl y™) 
AUPPM(KylyMI=ODPPM(KyigM) 
CONTINUE 

DO 4320 M=1,8 
ADPN(K,M)=UPNIK,M) 
ADPM(KyMQ=EPM(K,M) 
CUNTINUE 

CUNTINUE 

RETUKN 

END 


WRITE( 06,1313) 


(G{Ip),[P=1,4) 


{G(1IP),1P=1)4) 
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((ACTIP,JP),JP=1,4) ,IP=1,4) 


SUBRUUTINE NEWRAP(N,U,sDU,NPCT sOPSyMAXgMIgNI) 


IMPLICIT KEAL¥*¥8 (A-H,yV-Z) 


OIMENSTON UIN) MIC NPCT) PP 116)_GP(4),XP (4), AUX(16),LAUA(4), 
MAUX(4) ,ARLLodsuP(1lo),CP(16) UN) 


DG 10 L=1L,NPCT 
MI (11) =4 

NRHS=1 

DG LJU J=leMAX 
CALL ADULG(Nyl,NI}) 


CALL PUTRSINFUIT MI AP si PyC Psy PPyGPy_ XP yNRHSs AUX gLAUX»MAUX) 


CALL ERRPCTI(NPCT) 
Iédc=l 
DQ i¥9 I=i,NPCT 


LT1i1=4*( 1-1) 
READ(14*I22) 
be hee ek 

OU 20 1=i1,N 
U(1)=U(1)+bUt1) 
She Qs 

SZ=Ue 
Si=St+bDUu( 1) *0DU(1) 
$22 5240 L 1 *Ubi 
b=0S On bs bZo2a 


WRITEC6e 2222) 


{DUCT ITLIFKR) » KR=1, 4) 


WIE LL) ¢ Lb babah det 


Nor 


, | ( Sah ; 
be 


‘{eei= Shs (Pe. i= be tdi, Help bo 


(TAelMeXAMe293¢ T24 Vi pUUell etd Vat 
(\~tie He Al Be aa Fa 
st PTRUBS (OLD RUA CA) 9Xe (P1993 eC 01099 ¢ TORI OD | 
(AUG LOLP Ie (a1) 50 eee Tan! 
AS 


(AU AMe KU Ade KUA ce CHA: UX g Veg eV asda tan he be! } 24 Tas 


UT 3a ESA 


Sr (Aedelil EU) 
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2222 FORMAT(IH ,@E11.3) 


Lis=i1 
DO 31 J=1,NPUT 
TPi=44( 1-1) 


Weite CLOT Pr") (UC LTI+KR),KR=1,4) 
S40 (teat i 
Petite bet PS) GU To 200 
100 CONTINUE 
WRITE(6,1000) MAX,t,(U(1),1=1l»N) 
1000 FURMAT(LHO,* CONVERGENCE NOT ACHIEVED AFTER® ,14,iX,* CYCLES »/y5X,y 
L MPRESEN? TRELATIVeE ERROR 15°%)E10. 39/9 DXs"* PRESENT RESULT US*>57, 
2 (4E2027)) 
Cae sox Tt 
200 KETURN 
cNv 
SUBRUUTINE EKRPCTINPOT) 
IMPLICIT REAL*S8(A-H, O-Z) 
DIMENSION ABC (i6)_XP 04) 5X14) XN(4) gp AUX14) 
[1l2=1 
DO L0U T=1,NPUT 
0 pS ak WS J=154 
XP({J)=U. 
LQ AN{J)=0. 
Povtec tere wou TO 20 
1i4¢4=1-i 
READ(i4°114) (XP{K),K=1,4) 
Wea LZ’) t2) {ABC(K) ,K=1,16) 
1 WS acon § Wi eg 2s o 
gto ee XMXCABC XP XP 94 op 4eh 7 AUX) 
29 READ(L14'1) LAK(K),K=1,4) 
Repo tice lle) (AGC{K),K=1, 16) 
Tig=Illetl 
CALL XMXL(ABC eK yX9404 91, AUX) 
Ir (iseGeNPoT) GG Wi 59 
1i4=I[+1 
READ(L4*114) ( XN(K),K=1,5 4) 
Rasy bo Bp adie i Way) (AGC(K),K=L,16) 
Tig=lict+l 
CALL XMX( AGC yXNg XN 4949 he AUX) 
5U REAU(IL*I) (AUK K) »K=1,4) 
DO o8 K=194 
X(KJV=SXPCK)D EXTKVEXNIK) 
bY AUX{K)=AUA(K)—-X1{K) 
100 WRITE (SF LCCC) (X(K) pK=1 54) pL AULX(K) ,K=1 ,4) 
1000 FURMATCLH »Lore4tl 2o4y LUA, 4E12.4) 
rE TUKN 
END 
SUSKOUT INE” INV(AgN, Deh yh) 
LMPL Led REALS eo atAH-HyO-2) 
Ces? TS THe Pe sO eer SUBREUELINGE 


| : iA reapeah 
eXGeNg CASIO Xi ePl, * ADT “4 GaVSiH OA Ton | 
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| bose | 
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(S-DeH-A) sts848 oe 
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(Shi 


SEARCH FOR LARGESI 
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she) 
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INTERCHANGE ROWS 
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DAV lDEe (COLUMN BY MINUS PIVOT 
CGNTAITNED IN 


45 
46 
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46 


50 


DIMENSTON Alid,LC1),MC1L) 


D=Le 
NK=—-N 
DU 59% 
NK=NK+4N 
iL(K)=K 
M(K)=K 
KKAENK+K 
BIGA=A(KK) 
DO 20 J=KyN 
TZ=N*(J-1) 
DG 20 T=Ky,N 
IJ=1Z+1 


K=L_N 


TF (DABS EIGAI-CABS(A(TJ))) 


BIGA=A(1J) 
L(K) =i 
M(K)=J 
CONTINUE 


J=L{K) 
TF(J-K) 


> KI=K—-N 
DB K30 I=1,N 


KI=KI+4N 
HOLU=—-A(KT) 
J1l=K I-K+J 
A(KI)=aAlJT) 
A{ JI )=HCLOD 


I=M(K) 

Li feleK0 
JP=N¥*(1I-1) 
DO 40 

JK=ENK+S 
JI=JP+J 
HOLD=-A(JK) 
A( JK) =A(JT) 
A{JI)=HCLD 


J=1,1N 


ELEMENT 


35493 59 2D 


45945438 


BIGA) 


TF(VABS(EIGAI-1.D-20) 


D=0 ot 


WRITE{o,10GU) 


FOKMATCLHO, SUE TERMI NANT 


CALIi« tXelal 

DG 55 [=l,N 
leit ieKa 
IK=NKtI 


Be 5 aioe 


ALI K)=AC1KI/(-E ICA) 


460,746,488 


15420520 
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(VALUE OF PIVOT ELEMENT I5 


OS eOseCi ceed al 24 


West 
oe, oP ear «GS raei- 


“49.0 = Tay [Maat ad: Oh 


= i “?) a 
‘s! : eat 
corel avert 
: in a iS : - 
. “ - 


(Pu 


55 MGONTINUE 
REDUCE MATKIX 
DU 6) T=l1,WN 
IK=NK+] 
HOLL=AT1K) 
IJ=I-N 
DQ od J=1,N 
TJ=1J+h 
IFi{1-K) 6U_705,60 
O FM JAK) G2ytDy,E2 
2 KJ=1J—-1+K 


A(T J) =HOLL¥ AUK JI4A 15) 


605 CONTINUE 
Ol VpEs Rowe By PIVCT 
KJ=K-N 
DU 75 J=1,N 
KJ=HK SEN 
IF (J-K) 10, F570 
7Q0 A(KJI=A(KIID/BIGA 
75 CONTINUE 
PRODULT OF PIVUTS 
. D=C*81IGA 


REPLACE sPIWite Bie RECDPROOAL 


A(KK J=1./B81G4 
BOQ CUNTINUE 


FINAL RUW ANG CCLUMN INTEKCHANCE 


K=N 
109 K=K-1L 

IF(K) BDO 9b 50» 1:05 
1Q>5 [T=L() 


Lotie-Kg 120,120,106 


LO8 JW=N*{K-1L) 
JR=N*{(I-1) 
DO thebo J=1,N 
JK=Ju+J 
HULD=A(JKH) 
JIHSRK+/ 
A{( JK) =-A(JI) 

110 A(JI)=HOLD 

120 J=M{(K) 


IF(J-#) MOUs UC. L.3 


125 KI=K-N 
Diy ABO l=l 5 
KI=KI+tN 
HOLD=A(K I) 
JI=aKI-K+tJ 
ACK1)=-A( JI) 

130 AC JI)J=HOLD 
GU 4d tx 00 

150 (KE FURN 
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END 

SUBRUUT INE PUTRSINEQ SMI ,AP BP yCPy PPy GP, XPyNRHS gAUXsLAUX yg MAUX) 
IMPLICIT REAL*®& {A-H,O-Z) 
VBVIMENSTICN LAUX(1)+MAUX{ 1) 
DIMECNSTUN MitidsAP(ids BPLi de CP(L),PP(1),GP(1),XP(1) sAUXI1L) 
Tii=i 

[lé=1 

Télei 

[T2c=l 

N=M1(121) 

NPLi=M1I(T21+41) 

NN=N¥N 

NNPL=N*¥NP 1 

NN RHSEN*NRHD 

READ(L2* 112) {3P(K) sK=1_,NN) 
Pigedizei 

READ (IZITi2) (AP(K) »K=1L,NNPL) 
Il2Z=J12+1 

READ(GitH21) (GP(K) sK=i,yNNRHS) 
T1iL=ILiti 

CALi INV{ GP yNyCETyLAUXsMAUX) 
CALL XMX LBP, AP PP yNyNgNP1ly AUX) 
CALL XMX(UBPyGP yXPaNyNagNRHS y AUX) 
WRIZELGS* 121) (PP(X),K=i1,NAP1) 
I21=[21+1 

WRITEL L4¢1 22) (Ar (K) »K=L_9 NNRHS) 
I22= 122+1 

PEANEUSEO SZ3GC TO iol 
NEG@L=NEQ-1 

DU 100 T=d,NEQ1 

NMi=MI(1I-1) 

N=MI(1) 

NPL=MI(1I+1) 

NNML=N*NM1 

NN=N AN 

NNPL=N*NPL 

NN RH S=N*NRHS 

Sb AD hile ii 2) (CP(K) »K=i,yNAM1) 
fig=lLiatri 

READ Gir 8 Liz) (EP(K) sK=isgNnQ 


Iub2=Ub2+k 
READ(aZ'ITiz) (AP(K) ,K=igNNP1Q 
lizsbLtert 
READ(11' 111) (GP(K) sK=1ly,NNRHS) 
Lil=1il+l 


CALE XSXMX (Py CP »pPR yp bPaNoNM1l_9\yAUA) 
CALL INV GP yN,LET,LAUX,MAUX) 

CALL AMX (BP yAP yPPyNyNagNPL»_ AUX) 

Call ASXMX(GPelPsXPyGPeNeoNML yg NRHS AUX) 
CALL XMX (BP GP XP yNgNegNKHSs AUX) 
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WK EP ECL S Tie) (PP iK),K=i,NAPI) 
121i=121t1 

WRITE i4' 1e2) {XP(K) ,K=isNNRHS) 
I2é=122+1 

CONT TINUE 

NML=MI(NEQ-1) 

N=M1 (NEw) 

NNM i =NeNM 1 

NN=NAN 

NN RHS=NENKHS 

READE TZYI b2) (CP(K) sK=LyNNML) 
lizg=fig+t+l 

READ CiZ*hi2) (EPIK) sK=HLyNNI 
PY2='Pc+i 

READ TALT ILL) (GP{(KJ,K=i,NNRHS ) 
T1li=L1ii+l 

CALL XSXAX(GP CP yPP yg tPyiN,»NMLyNyAUX) 
CALL INV{(BPyNeLET,LAUX,MAUX) 

CALL XSXMX{(GP CP XP _oGPyNgNBMLygNKHS,ZAUX) 
CaLlL XMX(EP,GP yXP_yNyNygNKHS,AUX) 
WRITt(i4' 22) (xXP(K)_)K=1,NNKHS) 

DE* 290 T=deNE4 

T2i=NEQ-I41 

i22=I121 

N=MI(Te¢it+1) 

NM1=M1i¢t1I21) 

NNML=N*NML 

NNRHS=NMLENRHS 

REAVLListIi2i) {PP(K) ,K=iyNNML) 
READ(L¢' 122) (GE(K) »K=L»yNNRHS ) 
CALL XSXMK(GP»2 PP yg KP y XP yNML gp NagNRHS AUX) 
WRITE(14"1T22) (XP (KK) eK=1,NNRHS) 
CONTINUE 

RETURN 

END 

SUBRUUTINE KASKAPX (Ay by Capel »pMygh,AUX) 
IMPLICIT REAL*&S (A-H,O-Z) 

DIMENSIGN All) ,EC1L)-C01j,001),AUX(1) 


C MATRICES AKE GIVER GY ROWS 


100 


LN=L*N 

DO 200) 121,08 
NR=( 1-42) /N41 
NO=IT-{NKR-1) ¥N 
AUA(1) =U. 

Dis LJGO K=i,M 
Ni=M*(NK-1) 44K 
NZ=NC4+(K—-1 ) ¥N 
AUAL TJ =AUXTIT)+ECNLIFC(NE) 
CONTINUE 

DG 2090 T=1,LN 


CLMAK Lamy U8 
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DELP=ALTI~AUATT) 

CONTINUE 

RE TURN 

END : 

SUBROUTINE XMX(AsGeC ab pMyN yAUX) 
IMPLICIT REAL*¥S8 (A-H,0O—Z) 
DIMENSIGN AC1) 5201) 2C01),AUX01) 


C MATRICES ARE GIVEN BY ROWS 


100 


200 


LN=L*N 

DG 100 IT=1,LN 

NR= (1-1) /N+t1 
NC=I-(NR-i) *N 
AUX(1)=0. 

BO 100 K=i,M 
Ni=M*¥{NR-1) +K 
N2=NU+(K-1) *N 

AUX CT)=AUXTTIF4(INLIFBINZ 
CUNTINUE 
DO 2OU i=1,LN 
C(1)=AUXC1) 

CONTINUE 

RETURN 
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APPENDIX E 


CALCULATION OF NET CHORD ROTATIONS 

The top chord compression member of joists AX05 and AXO1 
have been investigated under various boundary conditions. Net chord 
rotations at both ends of the top chord compression member need to be 
calculated for the specification of displacement and spring boundary 
conditions. For displacement boundary conditions the net chord rotation 
is applied at both ends, while for spring boundary conditions the spring 
stiffness is calculated on the basis of the calculated net chord 
rotation. Calculation of the net chord rotations and spring stiffnesses 
are illustrated below. 
E-1 Calculation of Net Chord Rotation in Joist AX05 

The length of the top chord member 2 = 22.96" 

From Figure 5.5(b) and Table 5.1: 
. 5756908 


Vertical deflection at the right support (ug) 


Vertical deflection at the left support (u,) .4625646 


) 
) 

c) Relative deflection = .1131262 
) 


Rotation of chord due to this relative deflection = “ape 


.004927 
(e) Absolute rotation of the left support (8) = -.0109852 
(f) Absolute rotation of the right support (8, )= -.0021092 
(g) Net rotation from the chord at the left support (8';) 
= -.0109852 + .004927 = -.0060582 
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(h) Net rotation from the chord at the right support (6'-) 


= -.0021092 + .004927 = .0028178 


E-2 Calculation of Net Chord Rotation and Spring Stiffness in Joist AX01 
Length of the top chord member (2) = 23.17 inch 
From Figure 5.4(b) and Table 5.1: 

(a) Relative deflection = .5159207 - .4091547 = .106766 


(b) Chord rotation due to relative deflection gees = .0046 


(c) Net rotation from the chord at the left support 
= -.0110928 + .0046 = -.0064928 
(d) Net rotation from the chord at the right support 
= -.0017483 + .0046 = .0028517 
Once these net rotations are obtained,the spring stiffnesses are 
obtained as shown below. 


(e) Spring stiffness at left support (Ke) 


_ Moment at joint 5 for member 34 
Net rotation at the left support 


(f) Spring stiffness at right support (Ke) 


_ Moment at joint 6 for member 35 
Net rotation at the right support 
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From Table 5.2, 


_ 2205 7 “ey . 
Ke = ~Fo64028 * 339606 in-lb/radian 
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